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feat Capacity, Heats of Fusion and Vaporization, and 
Vapor Pressure of Tetrafluoroethylene’ 


George T. Furukawa, Robert E. McCoskey, and Martin L. Reilly 


The heat capacity of tetrafluoroethylene of 99.99)-mole-percent purity was measured 


in an adiabatic calorimeter from 16° 
5-degree intervals from 0° to 210 
fusion of three experiments gave 7,714.5 
142.00 0.01° K, 
temperature of 197.53 0.01 
the equation 


- 972.9810 
10€10Pmmtte = 4.71241 —— ms 


The entropy of the ideal gas tetrafluoroethylene at 
0.37 abs j deg~' mole~' or 64.55 
The entropy was computed from molecular and spectroscopic 


from the calorimetric data to be 270.06 
mole~'(1 cal= 4.1840 abs j). 


data to be 64.54 cal deg~' mole~'! 


l. Introduction 


In continuation of the program for the thermo- 
dynamic investigation of high temperature resistant 
polymers and their respective monomers or prototype 
molecules, heat-capacity, heats of fusion and vapor- 
ization, and vapor-pressure measurements were con- 
ducted with tetrafluoroethylene (C,F,). This sub- 
stance has recently become important as the mono- 
mer in the production of polytetrafluoroethylene 
(Teflon) plastic, which has exceedingly high chemical 
and thermal stability. Heat-capacity studies with 
this polymer have recently been reported from the 
National Bureau of Standards.’ 


2. Apparatus and Method 


The heat-capacity and heat-of-fusion measure- 
ments with tetrafluoroethylene were made in an 
adiabatic calorimeter which has been previously 
described.’ The heat of vaporization was determined 
in another adiabatic calorimeter constructed by R. B. 
Scott of the Bureau. This calorimeter was similar in 
design to those described by Osborne and Ginnings * 

nd by Aston, et al,® in which a throttle valve was 
ontained within the space enclosed by the adiabatic 
hield. Briefly, during the vaporization experiments 
known electric power was supplied continuously, 
and the vapor was removed isothermally from the 
‘alorimeter by controlling the throttle valve. The 
material removed from the calorimeter in a measured 
ime interval, was collected by condensation in a 
smal! high-pressure cylinder (150 ml) and weighed. 


' This paper is based on the work sponsored by the Ordnance Corps, U. 8. 
tment of the Army 
. Furukawa, R. E. MeCoskey, and G. J. King, J. Research N BS 49, 273 
952) R P2364 
'R. B. Scott, C. H. Meyers, R. D. Rands, Jr., F. G 
Bekkedah!, J. Research N BS 35, 39 (1945) RP 1661 
‘N Osborne and D. C. Ginnings, J. Research N BS 39, 453 (1947) RP1841 
aoe Aston, H. L. Fink, G. J. Janz, and K. E. Russell, J. Am. Chem. Soc, 73, 
) 


Brickwedde, and N 


to 210 
K was constructed from the data. 
+7 abs j mole~' at the triple-point temperature of 
The three heats of vaporization mé@asurements at the normal boiling 
K gave an average value of 16,821 
results of the vapor-pressure measurements from 142° to 208 


+ 4.816562 > 


K. A smoothed table of heat capacity at 


The average heat of 


10 abs j mole~! The 
K can be represented by 
10°47? +- 3.958793 = 10°77? 


2.427347 


10-°T 


197.53° K and | atm was computed 


0 09 cal deg . 


In the vapor-pressure measurements the same 
calorimeter used for the vaporization measurements 
served as a thermostated container. No provision 
was made to stir the sample. The pressures were 
read by means of a mercury manometer and a mirror- 
backed glass scale. The readings were converted to 
standard mm Hg (g= 980.665 cm/sec’, temperature 
0° C) on the basis that the local gravity is 980.076 
cm/sec’. 

The International Temperature Scale*® was used, 
and below 90° K, a provisional scale’? was used based 
on a set of platinum resistance thermometers cali- 
brated against a helium-gas thermometer. 
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3. Material 


The sample was originally received in a high-pres- 
sure cylinder (300 g. to 1.5 ft.*) and inhibited against 
polymerization and other reactions. To remove the 
inhibitor, the material was passed through a train of 
gas-washing bottles containing concentrated sulfuric 
acid, and to remove traces of permanent gases, it was 
later frozen and pumped three times. About 155 g 
of the purified material was collected. As this sub- 
stance is considered to be hazardous, particularly in 
the uninhibited state, it was kept at temperatures be- 
low that of dry-ice at all times, except when being 
weighed. 

The purity of the material was determined in the 
calorimeter before the thermal measurements from 
the equilibrium melting temperatures * on the bases 
that no solid solution was formed and that the liquid 
solution followed Raoult’s law. The equilibrium 
temperatures and the corresponding reciprocal of the 

| fraction melted, F, are given in table 1. 


*H. F. Stimson, J. Research N BS 42, 209 (1949) RP1962 

’ H.J. Hoge and F. G. Brickwedde, J. Research N BS 22, 351 (1989) R P1LISS 

*G. T. Furukawa, D. C. Ginnings, R. E. MeCoskey, and R. A. Nelson, J 
Research N BS 46, 195 (1951) RP2i91. 





Equilibrium me lting te m pe ratures of te trafluoroeth y- 
le ne 


Mole fraction impurity =0.0460 4 T, °K 


Reciproe= | of 
fraction melted 
vs 


142. O12 
142. 0018 
142. le 
* 142. (4S 


Triple-point temperature 
142.00+0.01° K 
Purity, 90.9» mole percent 


* Extrapolated 


The data were plotted and the purity was com- 
puted from the slope of the curve to be 99.99) mole 
percent, The intercept or the triple-point tempera- 
ture of pure tetrafluoroethylene was found to be 
142.00 +0.01° K. The results of the premelting 
heat-capacity measurements made during the course 
of this study indicated that the purity remained 
essentially the same. (The temperature in degrees 
Kelvin was obtained from the relation: °K =273.16° 

'C. Whenever the temperatures are expressed in 
the paper to the fourth decimal place, the last two 
figures are significant only in the measurement of 
small temperature differences.) 


4. Heat Capacity 


The heat-capacity measurements were made from 


16° to 210° K. In the computation to obtain the 
net heat capacities, the observed gross-heat capacities 
were first smoothed, then the smoothed tare heat 
capacities were subtracted at the corresponding 
equally spaced integral temperatures. Curvature 
corrections (see footnote 8) were applied to the gross 
and tare heat capacities wherever significant. The 
vapor-pressure data obtained in this work and the 
liquid densities reported by Ruff and Bretschneider * 
were used to make vapor saturation corrections." 

The results of the heat-capacity measurements 
gave a precision of +0.02 percent in the liquid and 
from +0.02 to +0.05 percent in the solid range, 
except for the lowest temperatures. There is, how- 
ever, a fairly large uncertainty in the vapor-satura- 
tion correction to the heat capacity at the higher 
temperatures because the amount of material in the 
filling tube is not accurately determinable. Con- 
sidering these factors, particularly the latter, and 
other possible sources of error, the uncertainty in the 
final heat-capacity values given in table 2 is believed 
to be +0.2 percent. In table 2 the values below 
15° K were obtained by the extrapolation (see foot- 
note 8) of a Debye function fitted to experimental 
values between 16° and 30° K. 


allgem. Chem, 290, 173 
RP 1483 


anorg 1433). 


(1046 


*O. Ruff and O. Bre tschneider, Z 
H. J. Hoge, J. Research N BS 3, 111 


Tasie 2. Molal heat capacity of tetrafuoroethylene at 
lemperatures 


Molecular weight = 100.02, ° K < 273.16 


aha j deq * mole 


0.0 

i] 
> i 
7.44 

7 

2 Liquid 
91 

65 
tS. a 

4 


* Cyacd. is the heat capacity of the condensed phase under its vapor pre 


k xtr ipolated 


5. Heat of Fusion 


The determination of the heat of fusion involy 
continuous introduction of electrical energy from 
temperature just below the triple point to just aboy 
it and the correction for the heat capacity of th 
sample and container. Also, a small premelting 
correction was applied. The results of the measure- 
ments are given in table 3. The values given in t! 
fifth column are the total heats required to melt t! 
156.2694 ¢ in the claorimeter. Taking into accow 
the precision obtained, and other possible sources : 
error, the uncertainty in the heat of fusion is cor 
sidered to be +7 abs j mole. 

TABLE 3 Volal heat of fusion of tetrafluoroethylen 


100.02, mass of sample =156.26004 ¢, triple-point ter 


ture =142.00° K, °K =273.16°+°C 


weight 


Molecular 


Heat ca- 
pacity 


rem perature Heat in- 
interval put 


ahs j abs j 
1, 326.9 12, 053. 7 
1, 707.4 12, O52. 6 
1, 457.7 12, 052. 7 


aba j 
13, 374. 1 
13, 755.1 
13, 504. 5 


137.7930 to 143.7828 
136.3972 to 144.1500 
137.4175 to 143.9057 


Mean heat of fusion 
Standard deviation 


6. Heat of Vaporization 


The heat-of-vaporization experiments were ma 
at the normal boiling point (197.53° K) as det 
mined by the vapor-pressure measurements. 1! 
experimentally observed quantity y, the energy 
input per mole of sample collected," is related to t! 
“HN, 8. Oshon (1930 


N. 8. Oshorne, BSJ. Research 4, 600 RP 168. 
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eat of vaporization, 1, by the expression 
Sie 

L T\ i?” 
\° is the molal volume of the liquid, 7, the 
te temperature, and P, the vapor pressure. 
<perimentally observed quantities, the value of 
mversion term (second term to the right of the 
ty sign), and the heat of vaporization are 
in table 4. The molal volume for the con- 
n term was obtained from the densities of 
luoroethylene reported by Ruff and Bret- 
der footnote 9). Considering the pre- 
and other sources of error, the uncertainty in 
eat of vaporization is believed to be +10 
mole The heat of vaporization, computed 
the Clapeyron equation by using the vapor- 
ire equation obtained in this work, the liquid 


(see 


ties reported by Ruff and Bretschneider (see 


Volal heat of raporization of tetrafluoroeth yle ne at 
197 53° K 


l-atm pressure 


Molecular weight = 100.02, °K = 273.1 


dp 
yt 
d7 


ahs } mole a 4 mole 
16, SN4 16, 81 
16, GOL r) 16, 830 
16, 800 819 


Mean 
Standard dev 


footnote 9), and the Berthelot equation of state, 
amounted to 16,850 abs j mole This is in fair 
agreement with the observed value. The constants 
for the Berthelot equation were reported by Renfrew 
and Lewis.” 


7. Vapor Pressure 


As there was no provision to stir the sample in the 
calorimeter, several series of vapor-pressure measure- 
ments were made from 142° to 208° K going down 
the temperature scale as well as up to check the 
temperature equilibrium. The agreement in the 
results of the measurements, as shown in table 5, 
indicates that equilibrium conditions existed. The 
vapor-pressure equation fitted to the experimental 
values is given by 
71241- 972.9810 

1 


log, Puomig= 4 + 4.816562 


10-?7—2 347 < 10> *7? +-3.958793 X 10°7T° 
In table 5 are given the deviations of the calculated 
values, based on this equation, from the observed 
values. Although the deviations from the equation 
vary somewhat systematically, no attempt was made 
to obtain a closer-fitting equation. 

The normal boiling point was taken to be 197.53 
+0.01° K, although the equation gives 197.525° K 
and the difference in the rounding can be measured 
by the thermometer and by the manometer. The 
dp/dT at the normal boiling point is about 41 mm 
Hg per degree. 


2M.M und EF. I 


Renfrew 


Vapor pressure of tetrafluoroethylene 
} } y 


K =27 


le 





8. Experimental Entropy 


The experimental data were used to compute the 
entropy of the ideal gas tetrafluoroethylene at 
197.53° K and 1 atm. In table 6 are summarized 
the entropy computations. The Berthelot equation 
of state and the critical constants reported by 
Renfrew and Lewis (see footnote 12) were used to 
compute the gas imperfection. The uncertainty of 
t0.37 abs j deg! mole~' for the entropy of ideal 
gas C.F, at 197.53° K and 1 atm was obtained by 
statistically combining the uncertainty in the entropy 
of liquid C,F, at 197.53° K and in the entropy of 
vaporization. The uncertainty in the gas imperfec- 
tion correction was considered negligible. 

In terms of the thermochemical calorie (1 cal 
4.1840 abs j) the calorimetric entropy of the ideal 
gas tetrafluoroethylene becomes 64.55 +0.09 cal 
deg™' mole ~'. 


TABLE 6 Summary of the Molal ¢ ntropy of tetrafluoroethyle ne 


Molecular weight = 100.02, °K = 273.16°+ °¢ | cal = 4.1840 abs j 


Syer (Debye 

Sie ce 1a. SOlid, numerical integration 
ASia.e fusion, 7714.5/142.00 

SSia’ we wr, liquid, numerical integration 
Entropy of the liquid at 197.53° K 

OSie vaporization 16821/197.53 

Siw gas imperfection 


Entropy, ideal gas at 197.53° K and | atm 


WasuHIneton, March 13, 1953. 


'9. Entropy From Molecular and Spe tro. 
scopic Data 


The entropy due to the external rotation was oy. 


puted from the molecular constants taken fron: ¢| 
electron diffraction results of Karle and Kark 
which the C—F and C=C bond distances and ¢} 
FCF angle were given as 1.313 +0.010 A, 

+ 0.035 A, and 114 +2°, respectively. 
of inertia calculated were /,=152.9 10-” g-om 
[,=254.7 10-* g-cm’, and /=407.6 10 
gm’. The fundamental constants used in {| 
computation were essentially those given in ¢ 
National Bureau of Standards Circular 461." 


The vibrational contribution to the entropy wa: 


computed from the recent frequency assignments 
D. E. Mann, N. Aquista, and E. K. Plyler at ¢| 
Bureau. 

The entropy computed from the molecular a 
spectroscopic data is summarized and compar 
with the calorimetric entropy in table 7. An 
cellent agreement is shown in the two results. 


Molal entropy of tetrafluoroeth ylene at 197.43 
from molecular and spectroscopic data 


TABLE 7. 


Molecular weight = 100.02, °K =273.16°+°C., 1 cal=4.1840 abs | 


cal deg 
‘, translation $7. 67 
¥, external rotation 
*, vibration 


S, total, ideal gas at 197.53° K and | atm 


S, calorimetric 


L. Karle and J. Karle, J. Chem 
}). Rossini, K. 8. Pitzer, W. J 
tt, M. G. Williams, and H. G 


Phys. 14, 963 (1950 
raylor, J. P. Ebert, J. E. Kilpatrick, ( 
Werner, NBS Circular 461 (Nover 
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{ydrothermal Preparation of Some Strontium Silicates 


Elmer T. Carlson and Lansing S. Wells 


A study 
thermal conditions. 
positions, Were prepared 


38rO0-28i0,4H,O, 28rO-S8i0,-H,O, and 


was made of the formation of hydrated strontium silicates under hydro- 
Seven hydrated compounds, believed to have the following com- 
SrO-28i0,-H,0, SrO-SiO,-H,O, 28rO-28i0,-3H,0, 38rO-28i0, 3H,0, 
35rOSi0,-2H,0. 


The anhydrous compounds 


SrO-SiO, and 28rO-SiOQ, were obtained by hydrothermal synthesis at temperatures as low 


as 124° and 142°. C., respectively. 
silicates possess hydraulic binding 


No evidence was found that the anhydrous strontium 
properties. 


X-ray powder diffraction patterns and 


some optical data are given for the various hydrates, and possible relationships with the 


calcium silicates are discussed. 


l. Introduction 


During the course of an investigation of the 
possibility of preparing compounds of the hydro- 
varnet type described by Flint, MeMurdie, and 
Wells [1],' but with strontium substituted for cal- 
cium, it became necessary to identify the strontium 
S.ilicates that were obtained in some of the experi- 
ments. In the absence of any specific data on the 
hydrated strontium silicates in the literature, a brief 
study was made of these compounds, and incident- 
ally of the anhydrous silicates as well. 

The anhydrous strontium silicates have been 
deseribed by Eskola [2], who also worked out the 
phase diagram for the system SrO-SiO,. The meta- 
and orthosilicates, SrSiO; and Sr,SiQ,, respectivel 
were the only compounds found in this system. The 
metasilicate is apparently hexagonal, but forms a 
solid-solution series with a-CaSiO, (pseudowollasto- 
nite) which is monoclinic, and may therefore be 
monoclinic also. An X-ray diffraction pattern for 
the orthosilicate Sr,SiO, has been published by 
S0O'Daniel and Tscheischwili [3]. Heats of formation 
spof both silicates from the oxides have been reported 
Bly Nacken [4]. 

Preparation of tristrontium silicate has recently 
been reported by Nurse [5], who also gives the X-ray 
(iffraction pattern. The compound was found not 
to be isomorphous with tricalcium silicate. 

The hydrated strontium silicates appear to have 
been studied very little. Jordis and Kantner [6], 
‘btained a precipitate approximating a 1:1 molar 
ratio by the reaction between Na,SiO, and a soluble 
trontium compound, and it was assumed that the 
oduct was a hydrated metasilicate. Zhuravlev 
7] reports that Sr,SiO, possesses binding properties 
ut gives no details as to the composition of the 
ivdrated product. 


Figures in brackets indicate the literature references at the end of this paper 


2. Anhydrous Strontium Silicates 


Strontium metasilicate and orthosilicate were pre- 
pared from stoichiometric mixtures of SrCQO, and 
silicic acid (reagent grade) by heating for an hour in 
an electric furnace at about 1,425° C. The composi- 
tion was checked by analysis and found to be in 
close agreement with that calculated, but the micro- 
scope showed a little inhomogeneity in the meta- 
silicate, and more in the orthosilicate. 

The metasilicate preparation was considerably 
improved in homogeneity by fusion in an oxygen 
blast. This also increased the average crystal size. 
Crystals of better shape, however, were produced 
by hydrothermal means, as described later, and 
these were used for determination of optical charac- 
teristics. The interference figure indicated that the 
crystal is biaxial, with a very small optic axial angle; 
the symmetry, therefore, is not hexagonal. The 
refractive indices were found to be a= 1.597, y= 1.637, 
in good agreement with those given by Eskola [2]. 
The 6-index is close to 1.597. 

An X-ray powder diffraction pattern? for a hydro- 
thermal preparation of SrO-SiO, is given in table 1. So 
far as known, no pattern for this compound has hitherto 
been published. The sample prepared by fusion of 
the mixed oxides gave a similar pattern except that 
many of the lines appeared to be double, triple, or 
quadruple, instead of forming well-defined single 
peaks. The significance of this difference is not 
known. It might conceivably indicate the existence 
of two forms having slightly different symmetry, but 
no evidence of any phase transformation was ob- 
served on differential thermal analysis of the prepae 
ration. This is in agreement with the findings of 
Eskola. 

The X-ray pattern resembles that of a-CaO-SiQ,, 
indicating a relationship between the two as sug- 


? X-ray diffraction patterns reported herein were made on an X-ray Geiger 
counter diffractometer. using Cu Ke; radiation 
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gested by Eskola. This is supported by the simi- 
larity in optical properties. 

No attempt was made to prepare anhydrous 
tristrontium silicate. 


3. Hydrated Strontium Silicates Prepared at 
Atmospheric Pressure 


3.1. 3SrO-2SiO2-3H20 


The compound 3SrO0-2Si0,-3H,0 was obtained by 
slowly adding sodium silicate solution (12 g 
Na,Si0,-9H,0 in a small amount of water) to a boiling 
solution of strontium hydroxide (30 g of the octa- | 


74 


hvdrate in about 800 ml of water). Boiling wa 
tinued for an hour, and the mixture was kept o 
steam bath overnight. The precipitate, appar 
gelatinous at first, became completely cryst; 
The ervstals were very small, acicular or lath-s} 
and slightly birefringent, with negative elong; 
From the analysis of two separate preparations 
molar ratio was calculated to be SrO:SiO,:H 
2.96:2:3.26 * for the one, 2.97:2:3.17 for the o 
From this the formula 3SrOQ0-2Si0,-3H.O is ass 
to be correct. X-ray diffraction data are give 
table | 

Larger crystals of the same compound were su! 
quently obtained hydrothermally. Some of 
are shown in figure 1. They were found 
biaxial, positive, with refractive indices: a=} 
B=1.602; y=1.614. The optic axial angle was ver 
small, 

Differential thermal analysis of this compoun 
shows a large endothermic effect at about 340° ( 
corresponding to dehydration. Two small ey 


thermic effects occur at 820° and 890° C, suggesting 
the possibility that the loss of water involves 4 


break~lown to the oxides, which subsequent 
recombine to form anhydrous silicates. 

During the preparation of this material, th 
reactants were in such proportion as to provide 
considerable excess of Sr(OH),, hence it may } 
assumed that no more basic hydrate may be forme 
by this method. An attempt to prepare a less basi 
compound by reducing the Sr(OH), concentration t 
the same molar strength as the Na.SiO; resulted in a 
amorphous product. 


3.2. 3SrO- 2SiO2-4H20 


By boiling a mixture of silica gel and strontiw 
hydroxide solution for half an hour, a different com 
pound was obtained. This consisted of distinctiy 
dendritic or feathershaped crystals with modera\ 
birefringence. From chemical analysis, the mola 
ratio SrO:SiO,:H,O was calculated to be 3.04 :2:4.0 
from which it was concluded that the compositio: 
was 35r0-2Si0,-4H,O. The differential heating cur 
on this material is similar to that of the 3:2: 
hvdrate except that the dehydration break is doubk 
with the first peak at 230° and a larger one at 300° ( 
indicating that the water is driven off in two stages 
The exothermic breaks occur at 800° and 890°. |i 
is concluded that the 3:2:4 hydrate first loses on 
molecule of water, forming the 3:2:3 hydrate, whic! 
is then decomposed on further heating. The 
conclusion is supported by results of hydrothermal 
experiments to be described later. 

The 3:2:4 compound was also obtained in a larg 
number of hydrothermal experiments. Figure 2 © 
a photomicrograph of one such preparation; th 
crystals are larger than those obtained by boiling 
but exhibit the same habit of growth. 


3 The percentage of loss in weight on ignition at 1,100° was taken as the perct™ 
age of HzO, with no correction for the slight amount of carbonate that may * 
assumed to have been formed from the CO; in the water and air in the react 
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ystals of 335rO0.28i10. 3H 2) 
with 38rO-SiO.-H.O 


Magnification, 1 


prismatic er ystals 


faint needles 


Comparison of the X-ray diffraction patterns of 
ver 20 preparations in which this compound was 
lentified showed that they were about equally di- 
vided into 2 types. The pattern given in table 1 is 
typical of the group arbitrarily designated type 1. 
In type 2 about half the lines were in substantially 
the same positions as in type 1, the rest being shifted 
distinetly to the right (that is, toward smaller spac- 
ings), except for the line at 2.07, which was shifted 
to 2.09. The significance of this is unknown. No 
correlation was observed between the X-ray pattern 
type and chemical analysis, presence of other phases, 
temperature of reaction, or refractive indices. 

The indices also varied among the several prepa- 
rations, but were difficult to determine precisely. The 
high index ranged from 1.627 to 1.637; the low index 
from 1.595 to 1.617. 

The same compound was obtained by heating 
siliea gel with strontium hydroxide solution on a 
steam bath overnight (temperature about 75° C), 
but the reaction was incomplete. Similar reactioa 
mixtures held at 30° and at 50° for several months 
failed to develop any of the 3:2:4 compound. On 
the other hand, the product was not entirely amor- 
phous in either case. An X-ray diffraction pattern 
of the material from the 30° test showed two rather 
diffuse lines. The 50° product showed the same 
lines somewhat sharper with the addition of a few 
others. Apparently conversion was only partial, 
and the erystalline product has not been identified. 


4. Hydrothermal Reactions 
4.1. Procedure 


Strontium hydroxide and silica gel were the raw 
materials for most of the hydrothermal experiments, 
lhe strontium hydroxide was “CP” grade and con- 
tained a variable amount of carbonate, which was 
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ystals of 35rO- 2810. 


HHO 


Fiaure 2 Dendritic er 


Magn 


fication, 700 


removed by a process of solution, filtration, and 
recrystallization. Four different lots of silica were 
used: (1) reagent grade “Silicic Acid’? containing 
0.16 percent of nonvolatile residue after treatment 
with hydrofluoric acid; (11) a commercial silica gel, 
extremely fine and highly reactive but containing 3 
percent of nonvolatile residue, including 1.17 percent 
of Na,O; (111) a commercial silica gel containing 0.45 
percent of nonvolatile residue; (IV) a specially puri- 
fied gel, relatively coarse and unreactive, containing 
0.10 pereent of nonvolatile residue. For other ex- 
periments, some of the strontium silicates were used 
as starting materials. The materials were mixed in 
the desired proportions, and a small quantity (usually 
about a gram), together with 1 to 2 ml H,O, was 
placed in a small platinum crucible supported in a 
pressure bomb containing sufficient water to main- 
tain a liquid phase at the desired temperature. The 
bomb was then kept in a controlled furnace for a 
period of several days. After the warming-up period, 
the maximum deviation from the reported tempera- 
ture was about +2° C. Upon removal from the 
furnace, the bomb was allowed to cool in front of a 
fan for about half an hour. Because of the slowness 
of the reactions, no significant changes are believed 
to have occurred during this cooling period. 


4.2. Results of Hydrothermal Experiments 
a. General 


Table 2 gives the results of a number of hydrothet- 
mal experiments, selected to present the essential 
information. Confirmatory or supporting data have 
been omitted in the interest of brevity. The experi- 
ments are tabulated in the order of increasing basicity 
of reacting mixture, and will be discussed in that 
order. There are certain observations, however, 
that apply to the study as a whole. 









Tapie 2. Results of hydrothermal experiments in the system SrO-SiO,-H,O 


Tem 
pera- Time Product 
ture 


Starting material 


Days 
SiO>-nHyO; 38rO0-28i07-4H20; SrO-SiO; 
SiO; nH,O; SrO-28i0, H,O 
Cristobalite; SrO-2810;-H;0 
Do 
Do 
Do 
38r0-28i0;-4H,0; SiO; nH,O 
SrO-28i10;- HO: cristobalite; SrO-SiOs. 
SrO-2810>-H,O; cristobalite 
Do 
S8rO.28i0;-H,0; cristobalite; SrO-SiO>. 
Do 
38r0-28i0»-4H,0; SiO;-nH,O 


som eSH,O+asio 
18r(OH )p 8HyO+S810, 
18r(OH), 8H,O0+4810, 


“inane 


18r(OH )p8H,:O+38i0, 


zySace 


jassr OH); 8H,0+5810, 


J 2 SHyO4+ 2810 





oO. 
S8rO 2810, H,O; SrO-SiO;; SIO; nH,O 
Do 
SrO.-2810;- H,O: S8rO-SiO:: cristobalite 
SrO.2810,-H,O; SrO-S8i0, 
8rO-28i10;-H,0; SrO-SiO:,; cristobalite 
Do 
SrO-S8iO:; quartz 
38r0-2810;-4H,0; SiO»-nH,O 


28r( OH), SH,»O +3810 SrO-28i02-H,0; SrO-SiO, 


7 
7 
7 
10 
i 
4 
6 
4 
5 
3 
6 
7 
5 





| do. \ 
38r0-2810,-4H,O; SiO;-nH,0. d 
0 
38r0-2810;-3H,O; 38rO-28i02-4H20; Si0y-nH,O ; : 
8 e 38r0-2810;-4H,0; Si0,nH,O te) | 
18r(OH)18H20+1810 13 SrO-Si0y-H:O; SrO-Si0» Bhi 
38r0.-28102-4H,0; SiO; nH,O. A ' 
38r0.-28102-4H,0; SiO; nH,0. sam 
SrO-SiO, H,O; SrO-Si0, 
SrO-SiOy- H,O; SrO-SiO»; 38rO-28i02-4HyO (trace ah 
38rO0.Si0;-2H,0; SiO;-nH,O hil 
S 12 38r0-2810,-4H,0; 38r0-28102-3H,O0; SiO, nH,O —y 
18r(OH)28H20+1810 6 38r0-28i0»-4H20; SiO»-nH,O me litte 
§ BrO SiO, . 
Do hes 
SrO SiO, Do 3 thei 
38rO-28i10)-4H,O J . 
0 t he 
( x 38r0-28102-4H,O; 38rO-28i10;-3H,0 
® 2810)-4H30 war! 
38r(OH )y8H,0+2810 ( J 3s 2810;-4H,0; 28r0- SiO, nose 
} 38r0-28102-4H,O; 38r0-28i02-3H,O; SrO-SiO,-H,O i 
7 SrO-SiO»; SrO SiO»-H,O; 28rO-SiO; \ 
6 SrO-SiO,; 28rO0-Si0, 7 t} 
4 5 Do orun 
2 38r0-2810y-4H,0 ture 
$ 38r0-28i10,-4H yO; 38rO0-28i0;-3H,O 4 
{SrO.2810)-4H,0 4 38r0-28i0;-3H,0 gene 
’ : 3 Do ‘ 
4 38r0-2810»-3H,0; 28r0-Si0, on 
4 38r0-28i10>,-3H,0; 28rO0-Si0,: SrO-Si0, ™ cate 
3 38r0-28i0)3H,0 
7 28rO0-2810,-3H,0; SrO-Si0, ugg 
ss 7 S8rO-SiO»-HyO; SrO-S8i0>» ; 
8SrO-2810;-3H30 6 SrO-Si0s; 38rO-28i0;-3H,0. _°! th 
7 28rO-SiO,; SrO-SiO; defir 
7 Do 
13 3Sr0-.28i0>-4H,0; Sr(OH);-8H,0 in | 
5 38r0-2810;-4H20; 38r0-28i0;-3H20; Sr(OH);8H,0 biax 
6 Do — 
& 38r0-28i0;-3H,0; 38r0-Si0»)-2H,O0: Sr(OH),8H,O m= COM| 
rer / 18 6 38r0-28i10;-4H,0; Sr(OH);8H,0 
28r( OH)» 8H,O+ 1810, 13 28rO-SiO»-HiO sugg 
& 38r0.28102-3H,0; 38r0-Si0,-2H,0; Sr(OH),-8H,0 Tl 
6 28r0-S8i0;-H,O0; 38rO0-28i0,-3H,O; 28rO0-SiO, 
6 25r0-SiO; narre 
uv Do 
28rO. 810, 7 Do and 
58r( OH )y 8H,O + 2810, 13 28rO0-SiO»-H,O; 38r0-28i0;-3H20; Sr(OH),-8H,0 : relat 
4 3S8r0-S8iO»)-2H,0; 3 SrO-28i0)-4H,0; 8r(OH)>8H,0 § 
7 28r0-8iO; H,O; 38rO0-S8i0;-2H,0; Sr(OH),8H,O 
: 4 38r0-8i0;2H,0; 28rO0-SiOsg, Sr(OH),-8H,O 
38r(OH), 8H,0+18i0; 7 38r0-SiO»)-2H,0; 28r0-8i0+~ 
4 38rO0.8i0;2H,0; 28r0-Si0, 
4 2S8rO-SiO,;; little unidentified material } | ' 
. 4 Do 
48r(OH);, 8H20+ 1810, 136 5 8rO-SiO,-2H,0; Sr(OH);8H,0 SrQ): 
ica 
Roman numerals indicate the particular Jot of silica gel used: (1) Reagent silicie acid, (Il) A commercial gel of high fineness and low purity, (III) A commercis ‘ } 
gel of medium fineness and purity, (IV) A rather coarse gel of high purity 7 anny 
com] 
ry 
he system was found to be rather complex. In | be that some of the compounds are metastable under By es 
addition to the 2 hydrated silicates described above, | all conditions, and have only a transitory existence. and 
5 , > > ys be . . 
5 new compounds were observed. All have been Another notable aspect of the system is that th poun 
tentatively assigned formulas, but the compositions | hydrated compounds occur only at relatively low fy“!!! 
unde 


| 
of some are not definitely established. This is due | temperatures. It is noteworthy that anhydrous 
to the failure to obtain pure preparations. It may | compounds were formed hydrothermally at tempers- 
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is low as 124°, and that most of the hydrates 
posed below 150°. Because of the narrow 
of stability, relatively small fluctuations in 
ature or in rate of heating may have affected 

sults 
vill be noted that the results are not always in 
|. Different products are sometimes obtained 
apparently identical conditions. It was 
ved that different types of silica may give 
‘ ent results. For example, the compound 
~-SiO,-2H,O appeared unexpectedly in a 2:1 
ive (experiment 64) in which the relatively 
se-grained silica gel 1V was used. In this case it 
be assumed that the low reactivity of the silica 
ted ina temporary excess of strontium hydroxide 
e the calculated ratio In other the 
rences are harder to explain. The type of silica 
| in each case is recorded in the table, although no 
mpt was made to study their differences svstem- 

lly 
\t temperatures slightly above 100°, the primary 
tion product in all appeared to be 
sr0-28i0,-4H,O, irrespective of the proportion in 
vhich the oxides were present. Crystals having the 
sume distinctive shape were frequently observed in 
action products that had been held at much 
higher temperatures, but it was found by X-ray 
liffraction and refractive index determinations that 
these had undergone decomposition without losing 
their over-all appearance. It must be assumed that 
the hydrate was formed initially, possibly during the 
warming-up period, and subsequently was decom- 
posed at the higher temperatures. 

\nhydrous strontium metasilicate (SrO-SiO,) and 
orthosilicate (2SrO0-SiO,) were formed at tempera- 
tures as low as 124° and 142°, respectively. In 
general, there was evidence of intermediate forma- 
tion of one or more of the hydrates, but the metasili- 
cate occasionally appeared in the form of spherical 
aggregates, as if it had formed directly at the surface 
of the silica grains. More often it occurred as well- 
defined tabular crystals of hexagonal shape, as shown 
in figure 3. These were found to be optically 
biaxial, with a very small optic axial angle. The 
compound thus is shown to be pseudohexagonal, as 
suggested by Eskola [2]. 

The rest of the products were obtained over a 
narrower range of oxide ratio in the reaction mixture, 
and will be considered in the following sections in 


cases 


cases 


s relation to that ratio. 


b. Mixtures Less Basic Than 1SrO: 1SiO, 


Hydrothermal treatment of mixtures ranging in 
SrO:Si0, ratio from 1:8 to 2:3 initially produced the 
3:2:4 hydrate. Continued treatment resulted in 
anhydrous SrO-SiOQ,, together with an unknown 
compound. In the highly siliceous mixtures the 
excess silica was converted to cristobalite above 180°, 
and to quartz in one case at 368°. The new com- 
pound was first detected by means of its X-ray 
diffraction pattern. It was subsequently observed 
under the microscope in the form of extremely fine, 
round, nearly isotropic grains, usually in relatively 


rystals of SrO-SiO formed under hydrotherma 


Ficure 3. ¢ 


conditions 


Magnificatior 


large aggregates. The mean index of refraction 
appeared to be about 1.574, but this figure is subject 
to considerable uncertainty Many of the grains 
seemed to have a much lower index, but this is 
believed to have been due to a shell of cristobalite 
surrounding a core of the silicate. The existence of 
solid solution appears unlikely in view of the con- 
stancy of the X-ray pattern. 

Because the new compound was not obtained pure, 
it was necessary to arrive at its composition in- 
directly. Mixtures of 1SrO:2SiO, yielded the new 
phase together with small amounts of both cristo- 
balite and SrO-SiO,. Increasing the strontium to a 
ratio 25rO:3Si0, resulted in the new phase plus 
SrO-SiO, but no cristobalite, indicating that the com- 
pound in question is richer in SiO, than the reac- 
tion mixture. For mixtures more siliceous than 
1SrO:2Si0,, the results are less easy to interpret. 
SrO-SiO, appeared in 1:3 mixtures, and even in one 
instance in a 1:8 mixture, at 150°. Bearing in 
mind the fact mentioned above, that the hydrate 
3Sr0-2Si0,-4H,O is formed initially even in the 
highly siliceous mixtures, it is possible that the sub- 
sequent decomposition of this hydrate leads most 
readily to the formation of the anhydrous metasili- 
cate. The latter is known to be relatively stable 
and, once formed, would have little tendency to 
react further; hence, it might occur in compositions 
far removed from a 1:1 ratio. Both cristobalite and 
SrO-SiO, were formed in small amounts from the 
2SrO0:5S8iO, mixtures, but the proportion of cristo- 
balite was estimated to be higher than that of the 
SrO-SiO,. Thus the ratio for the tiew compound 
appears to lie between 2:5 and 2:3; it is therefore 
assumed to be 1:2. Ignition loss determinations on 
several preparations indicated an average content of 
about 0.7 mole of water. Allowing for the anhydrous 
impurities present, the new compound is thus be- 
lieved to have the composition SrO-2SiO,-H,0. 
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c. Molar Ratio 1SrO: 1SiO 


Equimolar mixtures of strontium hydroxide and 
siliea gel, like the more highly siliceous mixtures, 
reacted to form 3SrO0-2Si0,-4H,O at the lowest tem- 
peratures used (110° to 114°). At 142°, and higher, 
the reaction product was always the anhydrous meta- 
silicate, SrO-SiO,. At intermediate temperatures, 
the nature of the product appeared to depend on the 
type of silica gel used, but the results were not always 
consistent even with a single type of gel. For 
example, table 2, experiments. 30, 35, and 36, reveals 
that both 3SrO0-2Si0,4H,O and 3S8r0-2S8iO)-3H,O 
were formed at 136° but only the former at 130° and 
137° in mixtures containing gel I]. The trihydrate 
was formed from gel ITV at 120° (experiment 27), 
but not from gel I at any temperature. The forma- 
tion of the highly basic compound, 3S5rO-Si0-2H,0, 
in experiment 34 has already been mentioned. 

In experiments 29, 32, and 33 a new phase appeared. 
This consisted of minute flakes with low birefringence 
and a mean index of 1.604. Repeated efforts to pro- 
duce this compound reasonably free from other phases 
met with failure. The best preparation, No. 33, con- 
tained 0.99 mole of H,O per mole of SrO-SiO,, while 
No. 32, with a higher proportion of anhydrous mate- 
rial, contained 0.74 mole of H,O. On the basis of 
this meager evidence the formula SrO-SiO,-H,O has 
been assigned to a compound. The X-ray pattern 
is given in table 

Treatment of alii SrO-SiO, at 132° left the 
material apparently unchanged, aah this is 


within the temperature range in which the various 
hydrates were formed from the hydrated oxides. 


d. Molar Ratio 3SrO:2SiO, 


Mixtures of strontium hydroxide and silica gel in 
3:2 ratio reacted as did the equimolar mixtures, ex- 
cept that the increase in SrO resulted in the formation 
of the anhydrous orthosilicate, 25rO0-SiO,, at the 
higher temperatures. This compound always ap- 
peared as irregular masses of small particles, often 
pseudomorphic atter 3SrO0-2SiO0,-4H,0. 
metasilicate, it was never found as well-defined crys- 
tals. In this, as in the preceding series, the results 
within a small temperature range are varied and in- 
consistent. 

In two series of experiments the starting materials 
were the 3:2:4 and 3:2:3 hydrates described above. 
From the results in table 2, it may be seen that the 
3:2:4 hydrate was converted to the trihydrate at 
temperatures of 132° C and higher. For reasons un- 
known, this result was not consistently obtained from 


mixtures of the hydrated oxides in the same molar | 


ratio. At 161° and above, anhydrous compounds 
were formed. 
The trihydrate behaved differently in that two 
ee hydrates were formed, at 124° and 138°, 
ectively. In the latter case, the compound was 
SrO.S SiO, H,O described in the previous section. 
The n new phase formed at 124° proved to be very 
elusive. It was obtained in small amounts in a few 
other experiments not included in the table, but 


Unlike the | 





never free of other phases. It appeared in the { pm 
of very fine needles with an average refractive in ley 
about 1.60. It could be distinguished from SrO-s.0 
H,O only by its X-ray diffraction pattern (see tsb), 
1). Although this compound was not obtained 
drothermally sufficiently pure for chemical anal 

& preparation obtained at a lower temperature, | 
method to be described later, was found to have 
following oxide ratio: SrO:SiO,:H,O=0.96:1:1 
The formula 2Sr0-2Si0,-3H,O is suggested as a | 
sibility on the basis of this analysis. 

It is worthy of note that although the 3:2:4 hy. 
drate was converted to the 3:2:3 by raising the tem- 
perature to 132° or higher, the reverse reaction was 
not obtained by holding the 3:2:3 hydrate at 106 
(experiment 55). 


e. Molar Ratio 2SrO: 1SiO. and 5SrO:2SiO, 


An increase in SrO to the 2:1 ratio gave rise to two 
new hydrated compounds. One of these, the tri- 
strontium compound, will be discussed in the next 
section. The other appeared, along with other crys- 
talline products, at 132° to 148° C. Analyses of 
preparations 66 and 72, which appeared to be most 
nearly homogeneous, gave the following molar ratios, 
respectively: SrO:SiO,:H,O=1.86:1:0.98 — and 
1.86:1:1.09. The formula 2SrO0-SiO,-H,O fits the 
data reasonably well. 

This compound occurred as small lath-shaped 
crystals, with parallel extinction and negative elonga- 
tion. It is biaxial positive in character, with low 
and high indices of 1.638 and 1.641, respectively 
The X-ray powder diffraction pattern is given in 
table 1. 

At 156° and above, the anhydrous orthosilicate, 
28rO-SiO,, was formed. Attempts to reverse the 
reaction by holding the anhydrous compound at 140° 
and below were unsuccessful. 


f. Molar Ratio 3SrO: 1SiO, and Higher 


Mixtures of strontium hydroxide and silica gel in 
3:1 ratio produced a new compound that has been 
referred to in the previous section. As in the cas: 
of the other hydrothermal preparations, it was al- 
ways contaminated with one or more other phases 
Analyses of three preparations gave the fo lowing 
SrO:SiO,:H,O ratios: No. 76, 2.9:1:2.2; No. 75, 
2.7:1:1.7; No. 80, 3.1:1:2.7. The first of these 
preparations contained a small amount of 2SrO0-Si0,, 
the second somewhat more. The third was prepared 
from an original 4:1 mixture and contained a small 
amount of Sr(OH),-8H,O, which probably accounts 
for the higher H,O content. On the basis of these 
analyses, it is concluded that the compound has the 
composition 3SrO-SiO,-2H,O. It was observed as 
very thin needles or plates (see fig. 1), often joined 
in spherulites. They were birefringent, with indices 
a=1.575, y=1.595, and showed parallel extinction 
and positive elongation. The X-ray diffraction pat- 
tern given in table 1 is a composite of five patterns 
on different preparations. 
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hydrate was formed together with 3S5rQ0-2Si0.- 
it 118° C. At higher temperatures the less 
ivdrate was changed to the anhydrous 2SrO. 
vhich frequently retained the peculiar shape of 
ginal crystals. The tristrontium compound 
med at temperatures as high as 183°, but at 
ind above only 2SrO-SiO, could be identified 
imple of the tristrontium silicate hydrate, after 
heated to 1,200°, was found to have decom- 
to 2SrO0-SiQ0, and SrO 


\ction of Water on the Strontium Silicates 


behavior of the various strontium silicates in 
resence of water was not studied thoroughly in 
nvestigation, but because of the possibility of 
iulie properties [7] it was considered desirable to 
a few exploratory tests. 
small quantity of the anhydrous metasilicate, 
SiO., was mixed with enough water to form a 

which was then packed into a vial and allowed 
io stand in a moist atmosphere for several days. No 
evidence of setting was observed. The orthosilicate, 
2Sr0-SiO,, behaved in the same manner, the only 
difference being that it became warm immediately on 
mixing with water. This may have been due to the 
presence of a little uncombined SrO in the prepara- 
tion used. This experiment fails to support the state- 
ment of Zhuravlev [7] that 2SrO-SiO, possesses hy- 
draulie properties. Negative results were likewise 
obtained when the pastes were kept at steam-bath 
temperature. Subsequent microscopic examination 
failed to reveal any change in appearance. As shown 


pin table 2, the anhydrous silicates also -were un- 


affected by water at temperatures just below those 
at which they are formed hydrothermally. 

The two anhydrous silicates were likewise unal- 
tered after being shaken with an excess of water and 
allowed to stand 3 days at room temperature. Titra- 
tion of the filtrates showed that a little strontium 


S livdroxide had gone into solution, but the concentra- 


tion was far less than the solubility of pure Sr(OH), 
at the same temperature. 

A single solubility test at 30° C was made on each 
of the two hydrated silicates that were produced by 
boiling; namely, 3S8r0-2SiO,-4H,O and 3Sr0-28i0,.- 
3H,O. The strontium hydroxide concentration was 
determined by titration after 1 day, and gravimetri- 
cally at the end of 10 weeks. There was no significant 
difference between the two values. The tetrahydrate 
residue was apparently unchanged during this 
period, and was in equilibrium with solution contain- 
ing 0.745 g of SrO and 0.304 g of SiO, per liter. For 
the purpose of comparison, the solubility of strontium 
hydroxide at the same temperature is about 10 g of 
SrO per liter. The trihydrate gave a solution of 

imost the same concentration, namely, 0.762 g of 
SrO and 0.303 g of SiO, per liter, but in this case the 
olid residue had been partially converted to the 
‘ompound tentatively designated 2SrO0-2SiO,-3H,O, 
nd the SrO:SiO, ratio in the residue had fallen to 
1.26. Another sample of the trihydrate, kept in con- 
act with water on the steam bath for 10 weeks, was 


’ 


completely converted to the same compound. The 
analysis of this material was the basis for the pro- 
posed formula. Progress of the transition was mani- 
fested by a curling of the acicular crystals, some of 
which developed a helical shape, as seen under the 
microscope. 

The glass containers used for the solubility experi- 
ments were strongly attacked by the solution, espe- 
cially at steam-bath temperature. This probably is 
the source of the additional silica required for the 
change from the 3:2:3 to the 2:2:3 compound. 


6. Relationship Between the Silicates of 
Strontium and Calcium 


It is of interest to consider possible relationships 
between the hvdated silicates of strontium and those 
of calcium. The latter have been studied extensively 
by a number of investigators, including Flint, Me- 
Murdie and Wells [8], and the subject has recently 
been reviewed thoroughly by Taylor and Bessey [9]. 

If the oxide formulas of the strontium silicate 
hvdrates described above are compared with a list 
of the reported calcium silicate hydrates, analogs 
will be noted for the I:1:1, 2:2:3, 3:2:3, 2::1, 
and 3:1:2 compounds. Whether any of these are 
actually isomorphous cannot be known without an 
analysis of the crystal structure. Visual comparison 
of the X-ray patterns of 3SrO0-2Si0,-3H,O — and 
3CaO-2Si0,-3H,O (afwillite) reveals considerable 
similarity. The refractive indices, birefringence, 
and crystal form also indicate the possibility of a 
relationship between these two compounds. To a 
lesser extent, resemblances may be observed between 
the compound designated 2SrO0-2Si0,-3H,O and the 
mineral crestmoreite, for which the formula 
2CaO-25i0,-3H,O has been given [8]. On the basis 
of refractive indices alone, it appears possible that 
the 1:1:1 and 3:1:2. strontium silicates may be 
related to the corresponding calcium silicates re- 
ported by Flint, MeMurdie, and Wells [8]. In the 
case of the 2:1:1 compounds, the relationship, if 
any, is probably with the calcium silicate having a 
refractive index of 1.64, sometimes designated the 
y-hvdrate. In these cases the data on crystal form 
are lacking, and the X-ray patterns are too complex 
for superficial comparison. 


7. Summary 


Hydrothermal treatment of mixtures of strontium 
hydroxide and silica gel in various proportions and 
at various temperatures produced seven hydrated 
strontium silicates, as well as the anhydrous meta- 


and orthosilicates. Table 3 gives the probable 
compositions of the hydrates, and their optical 
properties, so far as the latter could be determined. 
The X-ray powder diffraction patterns, together with 
that of SrO-SiO,, are given in table 1. Formulas 
given for the compounds having the oxide ratios 
| 3:2:3 and 3:2:4 are believed to be correct, and the 
| 2:1:1 and 3:1:2 formulas probably are correct, 
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TABLE 3 


Composition 


210 
2310 
2310,3H0 
2810.40 
10.H 


8S10,2Ho 


while the other 3 are somewhat in doubt. Both 
anhydrous compounds were formed hydrothermally 
at temperatures below 150 The hydrates 
38r0.2S8i0,.3H,0 and 38r0-28i0,4H,.0 were pre- 
pared at boiling temperature as well as under hydro- 
thermal conditions 

The anhydrous silicates did not appear to be 
strongly attacked by water, and no evidence of 
hvdraulic binding properties was found 


The X-ray patterns and photomicrographs were 
prepared by Barbara Sullivan and G. M. Ugrinic, 
thermal analyses were made by E. 5S. Newman, 
and many of the chemical analyses were performed 
by T.J.Chaconas. The writers are indebted to these 
people for their valuable assistance 


Character 


Some optical properties of the strontium silicate hydrates 
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Irregular 

Flat prisms 
Dendritic 
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Energy Reflected 


51.N 


From the Moon 
Ralph Stair and Russell Johnston 


tesults are given of some measurements on the ultraviolet and short-wavelength visible 
spectral radiant energy reflected from the surface of the full moon, made from October to 


Washington, D. C 
in direct 


December 1952, at 
the Fraunhofer bands as found 


Although the reflected lunar spectrum contains all 
sunlight 


with approximately the same relative 


intensities in the visible spectrum, intense absorption occurs for some of the ultraviolet 


wavelengths Selective absorption for wavel 


‘neths in the spectral regions of 380 to 


390 


millimierons and less than 360 millimicrons indicates the possibility of a lunar reflecting 


surface similar to that of powdered glassy sil 


1. Introduction 
\lan has long speculated about the moon 

n, its surface features, and its path in space. 

interest has stimulated the search for facts 
irding the exact nature and origin of the surface 
norama visible through the telescope 
however, generally agreed that the lunar 
features have been sculptured by cata- 
strophiec agents (either meteoric o1 voleanic, or both) 

to 10, inelusive].!. The resulting surface features 
liffer greatly from anything on the earth, except for e 
resemblance the few known terrestrial 
meteor craters. Because gravity on the moon 
only about one-sixth that on the earth and there 
s an absence of an atmosphere, the lunar craters 
are probably 25 to 50 times as large [11] as would 
result on the earth. 

\s the result of the absence of an atmosphere and 
moisture and, hence, of the usual types of weathering 
and erosion, the moon has retained records of many 
of its early catastrophic experiences. During its 
history about 16 times as many meteorites have 
collided with the earth, but their records have been 
large ly erased [1], unless we conclude that the 
encounters with the larger ones resulted in some of 
the geologic transitions indicated by abrupt changes 
between certain layers of the earth’s strata. 

The surface of the moon does change, however. 
It is affected by the sun’s rays, by gravity, and by 
tidal forces, by the temperature change from about 
F [15 to 18, inclusive] during the lunar day to 
about —150° F during its night, and by attrition due 
to falling meteorites, estimated at over one million 
per day [2]. These effects combine to produce a 
pulverization of the surface layer, which acts as the 
efficient insulating surface indicated by the character 
of the temperature changes on the moon's surface 
observed during solar eclipses [19]. 

The blackening of the old surface areas, or ‘“‘maria”’, 
may be due in part to exposure of the surface ma- 
terials to high frequency (short wavelength) radiant 
rgy. Oxygen, ozone, and other components of 
earth’s atmosphere act as a blanket that prevents 
earth from receiving ultraviolet and other solar 
ant energy of wavelengths shorter than about 
mu. Laboratory experiments show that many 
‘talline and glassy substances on the earth darken 

exposure to wavelengths within the spectral 
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icates 


range incident on the moon. The collection of a 
thin laver of fine meteoric material containing iron 
and other dark substances may be expected also to 
darken the lunar surface. 

The study of the moon through measurement of its 
effect upon reflected sunlight may be approached 
from several angles: through changes in intensity or 
polarization [11 to 14, inclusive] of the reflected 
radiant energy as a function of the angle of incidence, 
or through changes in the reflected spectrum caused 
by the lunar surface. This report deals with the 
integrated ultraviolet spectral intensities reflected 
for the fixed angular incidences corresponding to 
near full moon and with the moon near its most 
northern position in the sky. 

The observed relative spectral distribution of 
ultraviolet radiant energy reflected from the moon is 
very similar in quality to that emitted by the sun 
itself. All the Fraunhofer lines appear and with 
approximately the same relative intensities. Any 
differences in the two spectra result from selective 
optical absorption by the lunar surface. Except for 
a slight vellowing of the lunar image and for varia- 
tions in the ultraviolet and visible spectrum over its 
surface features [20], nothing has been reported in 
the available literature noting any other differences 
between the two spectra. 


2. Instruments and Procedure 

The apparatus employed in this investigation 
consists of a Carl Leiss double quartz-prism mirror 
spectrometer, using an RCA 1P28 photomultiplier 
as a detector. The light beam is modulated at 510 
c/s, and the output of the photomultiplier is fed into 
a tuned amplifier [21] and recorder (see fig. 1). A 
siderostat was employed for reflecting the beam of 
light from the moon into the spectroradiometer in a 
manner similar to that previously used with sun- 
light [22]. 

No condensing lens or mirror was employed, so 
that the resultant measurement was that for the 
integrated surface of the moon. The spectral- 
energy-response characteristic of the complete in- 
strument, including the siderostat. and photomulti- 
plier, was determimed by using a special tungsten- 
filament-in-quartz lamp [21, 23], together with a 
number of optical filters to reduce the lamp energy 
for the various parts of the spectrum (see fig. 2) to 
values approximating that of the moonbeam at the 
spectrometer slit. The radiant energy from the 





lamp was reflected into the spectrometer by the 
same siderostat mirrors, so that the spectral energy 
calibration for the moonbeam reduced toa simple com- 
parison of the recorder indications in the two cases. 

The high sensitivity of the detecting and recording 
equipment pe mitted the use of relatively narrow slit 
widths (spectral width approximately | mg at 310 mg 
and 2 to 3 my at longer wavelengths). These values 
are comparable to those employed in previous work 
with sunlight [24], so that the Fraunhofer structure 
of the measured radiant energy in the two cases is 
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quite similar (see fig. 5), although different spe. (ro. 
radiometers were employed. 

Measurements were made during four nights ar 
the ends of October and November 195 2, when th 
moon was near its full phase and also near its 1 
mum northern position, hence near its highest 
tude at the latitude of Washington, D. C. | 
best data were od during the night of Novem. 
ber 30—December when the moon was not on) 
nearest its full AS but was also at the highes; 
altitude for any night during the series of measure- 
ments. Also, the atmosphere was entirely free of 
clouds and showed least dust or haze scattering 0; 
this night. Data on the lunar altitude and air mass 
for the four nights were calculated in the usual 
manner by means of the celestial triangle throug) 
the use of the pertinent data published in thy 
American Ephemeris and Nautical Almanac for 1952 
for the solar and lunar positions. The resulting 
data are charted in figure 3. 

3. Spectral Radiant Energy Reflected From 


the Moon 

The spectral radiant energy reflected from th 
moon depends upon the optical and other physica 
characteristics of the lunar surface. Changes in the 
solar radiant energy are admittedly small [25] for 
all wavelengths penetrating the terrestrial atmos. 
phere. Variations at the moon as a function of tim: 
may be considered insignificant. In view of th 
fact that previous observers have found a marked 
variation in the light reflected from the moon as a 
function of the angles of incidence and _ reflection 
(11, 15], a similar behavior might be expected for th 
ultraviolet rays. However, at present the amount 
of this effect is unknown. Furthermore, variations 
in ultraviolet intensities over the lunar surface ar 
known to be appreciable [20]. Future studies for 
specific areas of the moon and at various angles of 
incidence and reflection should be interesting and 
informative. 

Terrestrial atmospheric absorption further puodi- 
fies the lunar reflected radiant energy. Mean spee- 
tral values (for ascending and descending moon) for 

lunar altitude of 65 degrees (air mass 1.10) an 
given in figure 4. When spectral radiant energy 
data are taken over a range of ascending (or descend- 
ing) positions of the moon and plotted logarithm 
cally [24] as a function of air mass and extrapolate 
to air mass equals 0, the intercepts represent th 
logarithms of the spectral intensities outside th 
terrestrial atmosphere. The data illustrated 
figure 5 were obtained in this manner. 

In order to illustrate better the similarity betwee: 
the spectral radiant energy reflected from the moo: 
as compared with that emitted by the sun, the data 
obtained for the sun at Climax, Colorado, in Septem- 
ber 1951, are taken from the previous publicatior 
[24] and reproduced in figure 5. Although thes 
data were obtained with different spectroradiometers 
at different times and places, the similarity in th 
two curves is striking. This is partly because th 
dispersions for the two instruments were not ap 
preciably different. However, slight differences re- 





RADIAN 


eur’ 
im ¢ 
that 
the 


con 
likel 
hye 
silic: 
Spor 
splu 
mea 
mitt 
teri 
abs« 
mat 


the 
Sica 
1 the 
| for 
MOs- 
time 

the 
rked 
as a 
‘Lor 
r the 
ount 
Ons 
y are 


: for 
al 


10d I- 


spec- 


erg\ 
end- 
hmi- 
ater 
th 
the 


' 


100! 
data 
tem- 
LL1O! 
hese 
‘ters 
the 
the 
ap- 


4 
¥ 
F 
' 
& 
4 
» 
f 
f 





WASHINGTON, D.C 1952 





370 390 “0 430 470 490 S10 $30 


WAVELENGTH , MILL IMICRONS 
S pe ctral distribution of the 
from the 


450 


radiant energy reflected 
moon 
in greater scattering of the data plotted in 

gure 6. 

(n inspection of the lunar relative to the solar 
wiant energy curve (fig. 5) discloses greater differ- 

es between the twe toward the shorter wave- 
lengths. A quantitative plot of this ratio (fig. 6) 
vives the relative spectral reflectivity of the lunar 
surface. Much of the seatter of the data, as in- 
dicated above, results from slight differences in the 
dispersions of the two spectroradiometers, inasmuch 
as a close inspection of the individual plotted points 
discloses that the higher values result from ratios be- 
tween peaks on the curves, whereas the lower values 
are associated with the Fraunhofer absorption bands. 
Three important characteristics of this reflectivity 
curve are worthy of note. First, the curve decreases 
n ordinate value with wavelength, thus indicating 
that the lunar surface has a lower reflectivity for 
the shorter wavelengths. Second, the band at 380 
to 390 mu indicates selective absorption of the lunar 
surface materials. Third, the sharp cutoff beginning 
at about 360 mu may be considered indicative of some 
special composition. 

lf the three special characteristics of the lunar 
reflectivity curve are considered in terms of possible 
materials present, and other known factors about 
the moon, such as its albedo, polarization, and heat 
conductivity are kept in mind, it appears not un- 
likely that a yellowish glass-like composition could 
be responsible for the observed phenomena. Certain 
silica glasses [26] have an ultraviolet cutoff corre- 
sponding closely with the observed curve. In a 
splintered or crushed form they would reflect a 
measurable amount of radiant energy after trans- 
mittanee through an appreciable thickness of ma- 
terial. A small iron content would result in selective 
absorption at 380 to 390 my and would give the 
material a slightly vellowish color. A pulverized 
rlassy silicate lunar surface would be highly insulat- 
and would produce characteristics compatible 
| temperature measurements obtained during 
pse and with lunar phase changes [15, 16, 19}. 
low average albedo [7, 27 to 30, inclusive], 
ut 7.3 percent, corresponds closely to the ex- 
ed reflectivity from glassy material. Although 
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measurements of the albedo of the moon are not 
precise, most of them fall below about 12 percent. 
Hence, the surface of the moon may be composed, 
at least in part, of powdered glassy silicates. In- 
cidentally, the high percentage of SiO, in the earth’s 
crust might suggest the possibility of a terrain 
similar to that of the moon had not air, water, 
erosion, etc., been present, 

The observations on polarization at the surface of 
the moon by Wright [12] indicate reflection by a fine 
texture and “point to pumiceous substances high in 
silica, to powders of transparent substances and to 
quartz porphyries and possibly to trachytes and 
granites as the materials we see at the moon’s sur- 
face.’ Similarly, the relative spectral reflectivity 
curve for the moon obtained in the present investiga- 
tion points to the possibility that the surface ma- 
terials are, at least in part, composed of powdered 
glassy silicates. Further refinements in the observa- 
tions of this interesting satellite are needed, however, 
before definitive conclusions can be drawn. 


4. Atmospheric transmittance and ozone 
The atmospheric transmission curve depicted in 
figure 7 is plotted in the usual way in terms of the 
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logarithm of the observed transmittances of unit 
atmosphere (at Washington) for the different wave- 
lengths as a function of the wavelengths. This, in 
turn, is expanded [31] according to the function 
(u-1)"A~* of the Rayleigh law of molecular seatter- 
ing, 
T 32r*(u—-1)°H log « 
oy / aN y4 ’ 
in which \ is the wavelength of the radiant energy, 
and uw is the index of refraction of the atmosphere 
Since, for the zenith position the atmospheric 
depth, /7, and the molecular density, .V, are constant, 
the resulting plot of the logarithm of the atmos- 
pheric transmittances becomes a straight line in 
those spectral regions wherein the Rayleigh law of 
pure molecular scattering is applicable. In as much 


as appreciable ozone absorption occurs omy at 
wavelengths shorter than about 330 my, the data 


herein recorded are inadequate for use in ozone 
determinations Between 300 and 330 My (fig. 4), 
the observed intensities were extremely low and the 
instrumental noise levels relatively high. With 
certain improvements in the equipment, it is hoped 
to reach sensitivities adequate for use of the appa- 
ratus in ozone determinations at night. As an 
alternative, the use of a condensing lens or mirror 
may be advantageous in supplying sufficient radiant 
energy from the moon for this purpose 


5. Summary and Conclusions 

This report presents the first observed ultraviolet 
photometric curve of moonlight from data obtained 
primarily on a single night, although measurements 
were made on three additional nights when the moon 
was near its full phase and at high altitudes. De- 
spite the fact that the measurements were made 
through the dense blanket of atmosphere over a sea- 
level station, interesting information was obtained 
having a bearing on the composition of the lunar 
surface. Further nighttime measurements at higher 
altitude stations with improved equipment should 
result in additional information not only on lunar 
reflectivity but also on ozone concentration. 

The extremely high sensitivity of the equipment 
lends its usefulness in other fields of research, in 
particular to stellar investigations. With telescopic 
magnification sufficient radiant energy from many 
stars, and also from small areas on the moon should | 
be available to permit precise determinations of | 


ultraviolet spectral intensities. Preliminary 
already made of weak fluorescent sources, refle: 
from dull surfaces, and of radiant energy from sy, 
sky areas (even during rainfall or after sundo wp 
indicate a wide range of possible application for th, 
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Phase Equilibria in the System MgO-TiO, 


L. W. Coughanour and V. A. DeProsse 


The phase relations in the 


svstem 


MegQO-TioO, 


were studied by means of solid-phase 


reactions between the components and by observing the fusion characteristics of specimens 


within the system. 


The existence of three binary compounds in the system was confirmed 


From the data obtained, an equilibrium diagram for the system is suggested 


l. Introduction 


study of phase-equilibrium relationships in the 

m MegO-ZrO,-TiO, has been initiated in con- 

ion with a fundamental study of ceramic 

ctries.! A survey of the literature pertinent to 
the problem revealed several discrepancies in the 
published data on the subordinate binary systems. 
lt was, therefore, considered essential that a prelim- 
nary study of the binary systems be made. 

(n equilibrium diagram for the MgO-TiO, system 
was first published in 1932 by von Wartenberg and 
Prophet (1].° Two compounds, 2MgO-TiO, and 
\VicO-2TiO,, were postulated. The subsequent work 

other investigators [2] showed, however, that 
there are three binary compounds in the system. 
They are the two listed by von Wartenberg plus 
\icO-TiO,. This compound occurs in nature as the 
mineral geikielite. 

The existence of three compounds in the system 
was confirmed in the present work by means of solid- 
phase reactions between the components MgO and 
TiO... These reactions took place at various temper- 
atures between 1,200° and 1,550° C. The complete- 
of reaction and compound formation was 
observed petrographically and by means of X-ray 
diffraction patterns of the fired samples. 

The determination of the fusion point of each of 
the three compounds and of the solidus and liquidus 
temperatures at various points across the system 
supplied data from which a revised equilibrium 
diagram was constructed 


2. Sample Preparation and Test Methods 


The following materials were used as components 
n the preparation of samples for this study. 

MgO—Water-clear crystals of Norton Co. artificial 
periclase of nominal purity over 99 percent were 
selected. The crystals were ground with a mullite 
mortar and pestle to pass a 325 mesh screen. A 
spectrogram of the ground material showed only 
weak lines for the principal impurities, Al, Si, and T1. 

TiO,—Highly purified titania of nominal purity 
99.9 percent was used. This material was 
secured from the Titanium Division of the National 
Leal Co. A spectrogram indicated that the prin- 
cipal impurities, Si, Cu, and Al, were present in 
amounts less than 0.01 percent. 
Calculations of weight composition were made to 
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+0.01 percent, no correction being made for the 
percentage of purity of the components 

The components in sufficient amounts to makea 
10 g sample were weighed to the nearest milligram 
and mixed in the dry condition by shaking in a large 
container. They were then mixed with approxi- 
mately 5 percent by weight of a 5-percent soluble- 
starch solution binder, and 1-in. diameter disks were 
formed in a hardened steel mold at a pressure of 5,000 
lb/in.*. The pressed disks were calcined for 4 to 16 
hr at 1,200° C. During heating the disks rested on 
platinum foil. Oxidizing conditions were maintained 
by the use of platinum-wound or Globar resistor air- 
atmosphere furnaces. After cooling, the disks were 
pulverized using a mullite mortar and pestle, the 
powdered material was remixed with starch solution 
binder, and repressed at 5,000 lb/in’. The repressed 
disks were then refired for 4 to 5 hr at 1,350° C, again 
under oxidizing conditions. Petrographic examina- 
tion, supplemented by X-ray analysis, showed that 
solid-phase reactions between the components were 
essentially complete after these two heat treatments 

Following the second calcination the disks were 
reground, remixed with binder, and disks about ', in 
high were formed at 15,000 Ib/in2 in a ‘4s-in.-diameter 
mold. From these disks the actual fusion test speci- 
mens were ground in the form of small four-sided 
pyramids, grooved on each side [3]. The fusion tests 
and softening-range tests were carried out in a thoria 
resistor furnace [4]. The specimens were placed in 
the furnace on a disk of a platinum-rhodium alloy 
The metal disk was, in turn, on a support of ThO, or 
BeO. The heating rate during the tests was about 
3 deg C/min until the solidus temperature was ap- 
proached, when it was reduced to about 2 deg C/min. 
Oxidizing conditions were maintained by slight con- 
vection currents of air through the furnace. The 
temperature and the fusion characteristics of the 
specimens were observed with an optical pyrometer 
calibrated at the National Bureau of Standards 
The melting point of pure platinum wire was observed 
as a spot check on the pyrometer calibration. Agree- 
ment within 2 deg C of the accepted value of 1,769° C 
on the 1948 International Temperature Scale was 
obtained. 

In the determination of the solidus and liquidus 
temperatures of the test specimens, the first sign of 
liquid formation, rounding of the corners of the test 
pyramid, was recorded as the solidus temperature. 
The temperature of complete fusion, as observed 
through the optical pyrometer, was considered to be 
the liquidus temperature. 





In such determinations, there are several possible 
sources of error. These include the slight deviation 
from black-body conditions, the introduction of mi- 
nute quantities of impurities in the forming and grind- 
ing of the specimens, the slight reduction of the TiO, 
on heating, and the inherent difficulty in the visual 
observation of fusion characteristics 

In the oxide-resistor furnaces used for the fusion 
tests, black-body conditions were rather closely ap- 
proached. In order to render the grooved pyramid 
specimen visible in outline, it was necessary to have 
a slightly concave surface on the metal disk on which 
the sample rested during heating. This concave 
metal surface reflected less light to the pyrometer 
than the sample did, giving the disk the appearance 
of being slightly cooler than the specimen 

Partial reduction of the titanates was observed in 
compositions richer in TiO, than MgO-TiO, when 
the specimens were quenched or cooled rapidly. It 
is probable, therefore, that at high temperatures and 
during the fusion tests there was some loss of oxy- 
gen from the samples. This loss might be thought 
of as introducing a minute quantity of a third com- 
ponent, or impurity, into the samples. It is be- 
lieved that the error introduced into the fusion-point 
measurements by this reduction is not significant. 
Neither by petrographic nor X-ray methods was it 
r0ssible to detect any evidence of the formation of a 
— oxide of titania, even in specimens that had 
been quenched and were a dark-gray color. This 
would tend to indicate that the amount of redue- 
tion is small. Agamawi and White [5] conducted a 
weight-loss test on rutile at various temperatures up 
to 1,550° C. They reached the conclusion that the 
oxygen loss, on heating rutile in air, is practically 
negligible Likewise, Moore [6] showed that, for 


TABLE 1 I 


single erystals of rutile, less than 0.1-percent 
tion caused a marked color change. 


Probably the largest source of error in the f) sjop. 
point measurements was in the problem of sele +tin, 


the solidus and liquidus points by visual obs 

tion. The accuracy would depend upon the vis: 

of the melts, the composition, and its relation 
eutectic point or compound, and the experienc: 
judgment of the observer. In this particular sy 

the fusion of compositions richer than 2MgO.T 

MgO was difficult to observe due to the high vise 
of the melts. 


Considering these various factors, it is believe 


} 


that the recorded temperatures are correct to 


deg C. 


After the fused specimens had cooled in the fu. 


nace, they were analyzed with the petrographj 
microscope. Such examination 


various compounds in the system, and by the fa 


that large well-formed crystals were obtained. Thy 
optical properties of the compounds may be taby 


lated as follows: 


Approximate indices of 


Optical 
refraction 


Con ul 
ompound character 


1.736 
1.959 
1.95, w 


MgO Isotropic 
2MgO-TiO do 
MgO0-TiO Uniaxial 
negative 
Biaxial 
negative 
niaxial 
positive 


2.28 
MgO0-2TiO 2. 


TiO UC 


sston characteristics of compositions in the system Mat -TIO 


Effect of heating 


Fusion 
st irted 


Results of petrographic examination 
(Specimens examined were not 
quenched and the phases observed 
may not be those in equilibrium at 
the liquidus 


Fusion 
completed 


Mge0O-+ 
Do 
Do 
Do 
Do 


2Me0.TioO 


Do 
2Mg0.TiO;:+ MeO. TiO 
2M¢e0.TiO;.+MeO. Tid 
Do 
Do 


tract 


Do 
Mg0O-TiOe+ MeO-2TiOo(trace 
MgO. TiOe+Me0.2TiO 

Do 

Ik 


Do 
MgO. 2Ti0, 
MgO0.2TiO;+ TiO 
Do 
Do 


Do 

Do 

Do 

0 
MgO-2Ti0;.+ Tio 


was facilitated }y 
the wide difference in the optical properties of th 


flaa a anc itiniiins statin DdSiadin, ahead, indi 





Serta? fey ees 


the biaxial compound MgQ-2TiO., 2V was 
ted, from interference figures, to be 70° to 
Calculation of 2V from the indices of refrac- 
ive a Value of 68°. 

results of the fusion tests and of the petro- 
« analysis of the nonquenched specimens are 
n table 1. 


tesults and Comparisons With Previous 
Work 


3 


equilibrium diagram for the system as deter- 
| by von Wartenberg is shown in figure 1. 
constructed from the results of the present 
is shown in figure 2. Comparison of the two 
ams shows that, apart from the introduction 
' third compound, the major difference between 

e two is the lower fusion points obtained in the 

esent study. It is believed that these temperature 

rences are explained simply by the greater re- 
nent of the methods used in the study covered 
by this report. 

Other fusion-point data published on the system 

extremely limited in extent. The fusion point 
of the compound MgO0.2TiO, was determined by 
Sigurdson and Cole [2,e] to be 1,645° C under oxi- 
ing conditions. This is in excellent agreement 
vith the presently determined value of 1,652° C 
St. Pierre [7] determined the fusion point of highly 
uirified TiO, to be 1,840° 10° C, again under 
oxidizing conditions. This value, also, is in excellent 
wreement with the value 1,839° C observed in the 
resent study. 

Biissem, et al. {2,a] conducted an X-ray investiga- 
tion of the system They did not determine an 
equilibrium diagram, but they did obtain results 
that indicated, to them, the existence of a high- 
temperature solid solution of MgO in 2MgO-TiQ,. 


They selected compositions in the area between 
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MgO and 2MgO.-TiO, and heated them to approxi- 
mately 1,460° C. They then compared the X-ray 
powder patterns of the cooled samples with the 
patterns of mechanical mixtures of the same over-all 
compositions, the mechanical mixtures being pre- 
pared from the components MgO and 2MgQO.-TiQ,,. 
They noted that the diffraction lines of MgO were 
sharper and clearer for the mixtures than for the 
heated compositions. This effect they attributed 
to a high-temperature solid solution of MgO in 
2MgO-TiO, with an unmixing or exsolution reaction 
on cooling. The unmixing was thought to result in 
extreme subdivision of the MgO grains with the sub- 
sequent lack of sharpness in the MgO diffraction 
lines. 

A similar test was made during the present study 
on the composition 4MgO:TiO,. The results ob- 
tained were negative, there being no significant 
differences in the diffraction patterns of the heated 
sample and of the mechanical mixture. 

Several attempts were made in the present study 
to detect solid-solution formation in the system, 
Particular attention was given to the high MgO por- 
tion of the system, in order to verify or disprove 
the existence of the high-temperature solid solution 
of MgO in 2MgO-TiO, postulated by Biissem. In one 
method, samples high in MgO content (9MgO:TiOg, 
4MgO:TiO,, 3MgO:TiO,) which had been previ- 
ously calcined were used. Portions of the calcines 
were heated for 2 hr at 1,400° C and for ') hr at 
1,650° C and were cooled slowly. The cooled samples 
were examined petrographically and by means of 
X-rays. The petrographic examination showed that 
the index of refraction of MgO, and of 2MgO-TiO,, 
had not deviated from its normal value in either 
ease. In addition, the MgO grains did not have the 
appearance of having undergone an unmixing re- 
action. X-ray diffraction patterns showed no sig- 
nificant intensity variations or parameter changes. 


TEMPERATURE , % 


o 
> 


MQ MOLE , PER: 


lo the MeOTi-O system, 


Figure 2. Proposed diagram 





In another attempt to verify any solid-solution 
formation, specimens of the same compositions were 
heated to 1,460° or 1,500° C and the temperature 
was maintained constant for 3 hr. The specimens 
were then quenched in water. Again, petrographic 
examination of the quenched specimens showed no 
change in the index of refraction of either MgO 
or 2MgO-TiO,. The results of these tests indicate 
that there is no solid-solution formation between 
2MgO-TiO, and MgO. Index of refraction meas- 
urements of these two compounds were accurate to 
about 2 units in the third decimal (+0.002). With 
this degree of accuracy, solid solution formation to 
the extent of 1 percent could have been detected. 

Similar tests at various points on the diagram 
failed to indicate any solid solutions in the system. 
Finally, the fact that the solidus lines obtained are 
horizontal is indicative of the absence of solid-solution 
formation. 


4. Summary and Conclusions 


Twenty-six compositions in the system MgO-TiO, 
were prepared. They were studied from the stand- 
points of solid-phase reactions and fusion behavior. 
The existence of three binary compounds in the 
system was confirmed. They are 2MgQO-TiO,, 
MgO-TiO,, and MgO-2TiO,. Their melting points 
were determined to be 1,732°, 1,630°, and 1,652° C, 
respectively. Four eutectic pomts were located as 
follows: 21l-mole percent of TiO,, and 1,707° C; 
44-mole percent of TiO), and 1,583° C; 56-mole 
percent of TiO, and 1,592° C; and 86-mole percent 
of Tit » and 1,606° C, No evidence of solid-solution 









formation was obtained. From these dat 


equilibrium diagram for the system is suggest; 


This work was sponsored by the Office of Ord» ay, 
Research, Department of the Army. 
All X-ray analyses for the project were perform, 


by R.S. Roth. 
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Pairs of Normal Matrices With Property L 
Helmut Wielandt ' 


A short proof is given, under weaker assumptions, of the following theorem first proved 


If the eigenvalues a,;, 8 


that the eigenvalues 4 


bv N W legmannh 
bered in such a way 


for i l, 


n 


> 


to have property 
8.(k=1, may be 
d in such a way that the eigenvalues of 
8B are given by aa,+ 88, for all complex num- 
x, 8. Though every pair of commuting matrices 
property L, the converse is not true.2 However, 
{ and B have property Z and are normal, then 
BA according to a theorem of Wiegmann.’ We 
to show that the assumptions of Wiegmann’s 
We begin 


some 


nH n matrices A, B are said 


heir eigenvalues ay, n 


vem may be weakened considerably. 
roving the following lemma which 
est in itself, 
MMA Let A and B be arbitrary n 
the set Zz of all points in the com ple r plane for 
A ~B is normal. is either the whole plane, a 
right line. a circle, or it contains, at most, two points. 
lf B is normal, neither the circle case nor the two-point 
IfZ is the entire plane then AB BA. 


has 


n matrices. 


Ther 


} 


may enter. 


Proor: A+B is normal if, and only if, 


0). (1) 


‘ 
‘ 


J iy, (1) is equivalent to a set of 2n* real 
inear equations in z,y, and (2°+y' Each of these 

equations defines either the whole plane, a 
circle, a straight line, or else, at most, one point. 
The intersection of these 2n’ point sets obviously is 
one of the six types described in the lemma. That 
all these types may indeed occur is shown by the 
following matrices A+ 2B 


) 


5 
é 


B: 


0 


-_ 
In case B is normal, the terms of (1) containing 
cancel, hence Z is determined by linear equations 
inzsandyonly. So Z is neither a circle nor a pair of 
poin ts 
can University and University of Tibingen. This paper was prepared 
« National Bureau of Standards contract with The American University 
S. Motzkin and O. Taussky, Pairs of matrices with property L, Trans 
th. Soe. 73, 108-114 (19452 
Viegmann, Pairs of normal matrices with property L, Proc. Am. Math. 
6 (1953), 


of two normal n 


and all complex values of 2), 


2) 


B may 
by ¥ 


n matrices A, 
zh are given 


BA, 


nuh 


Y 


of C a- 
then AB 


In case Z is the entire plane (1) holds for every 
hence the coefficient of vanishes, that 
BA*—A*B=0. Since B commutes with A*, it 
commutes with every polynomial f(A*), hence, with 
A itself. (If A is normal, then A=f(A*) for some 
polynomial f. This well-known fact follows by 
transforming A into diagonal form by means of a 
unitary transformation.) The lemma being proved, 
we turn to the generalization of Wiegmann’s theorem 

Tueorem: /f the elge nvalues a,, By, of two normal 
matrices A=(a,,), B 
way that the ergenvalues 
ine qualitie S 


Is, 


also 


(h ) may he ordered in such a 
; A B satisly the 


Yi ol 


for some values » (A=], 1) which are the 
vertices of a polygon containing O in its interior, then 
AB=BA. (This obviously contains Wiegmann’s 
theorem since if A and B have property ZL then (2 
holds even with equality sign for every 2.) 

Proor: According to a theorem of Schur,! 
have 


} 


we 


i 


‘ 
; > 
i 


20 4) 


+ ‘ 
oe 


>, Where the equality holds if, and only if, 
From (3) and (2) we have for 


7) - 


for every 
A+2B is normal. 


~A 


ob ,, (4) 


aT 
P. a 


which reduces to 


“ 
S (a; 
‘ 


2p )>0 


s (a, ra 


i 


k zh ~h } 


k 


since Say ?=>) aq? and 5) by)?=35)8,? in view 
of the fact that A and B are normal. Since the 
left-hand side of (5) is linear and homogeneous in 
xz and y, (5) defines either a half plane containing 0 
on its boundary or the left-hand side of (5) vanishes 
identically. The former case cannot occur, for the 
half plane contains 2),...,2,, hence it contains 0 in its 
interior. So in (5), and hence in (4), equality holds 
for every Combining this with (2) and (3) we 


*1. Schur, Uber die charakterischen Wurzeln einer linearen Substitution mit 
einer Anwendung auf die Theorie der Integralgleichungen, Math. Ann. 66, 
488-510 (1909); Theorem II. The theorem of Schur has recently been used in a 
related problem by H. Schneider, Theorems on normal matrices, Quart. J. Math 
Oxford [2] 3, 241-249 (1952 





conclude that in (3) equality holds for z= 2, hence 
A+z2,B is normal. Since 2,....2; are not collinear 
and B is normal, the lemma asserts that A+ 23 is 
normal for every z, hence AB= BA. 

Remark: The theorem is best possible in the 
sense that there exist noncommuting normal matrices 
A, B such that inequality (2) holds for every 2 of a 
closed half plane containing 0 on its boundary (and 
for a fixed suitable ordering of the eigenvalues 
a8). An example is given by 














with the ordering a,=1, 
Here we have y,(2)=1, 
hence >})\y,(2)/?2> Sola 
(though y,(2)=a,.+ 2B,, 
only). 
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Heat Capacity of Gaseous Hexafluoroethane’ 
John S. Wicklund, Howard W. Flieger, Jr., and Joseph F. Masi 


An accurate flow calorimeter has been used to measure the heat capacity (C, 
+ 10°, 
The results are believed to be accurate to 


hexafluoroethane (C,Fs) at 50°, —20", 


1.5-atmosphere pressure 


of gaseous 


+ 90° C and at 0.5-, 1.0-, and 


0.1 percent 


+ 50°, and 


The values of C, have been extrapolated to zero pressure at each temperature; the ideal-gas 
heat capacities thus obtained at the five temperatures are, respectively, 20.99, 22.79, 24.43, 


26.45, and 28.24 calories mole~! degree™'. 


Calculated results, using recent frequency assign- 


ments and molecular data, were 0.4 to 0.7 percent higher than the experimental results 
The values of the pressure coefficient of heat capacity at the five temperatures have been 
used in conjunction with literature data to obtain an equation of state for hexafluoroethane 


at low pressures 
l. Introduction 


irate measurements of gas heat capacity over 
- of low pressures and at several temperatures 
‘important information about both the ideal 
nl the low-pressure equation of state. The 
alorimeter used in this investigation has been 
Spreviously deseribed.2 * This measures the heat 
capacity of gases with an accuracy of the order of 
10.1 percent, from —50° to + 100° C, and up to about 
> atm 
© The heat capacity of gaseous hexafluoroethane is 
Sof particular interest because of the question of the 
Ssize of the potential barrier to internal rotation. 
Pace and Aston * made calorimetric measurements 
Bon the solid and liquid, and obtained third-law 
entropies of the ideal gas; comparison with certain 
pectroscopic calculations led them to a value of 
50 calories mole! for the potential barrier. 


flow 


2. Experimental Procedure 


2.1. Material 


A specially purified sample of C.F, was obtained 
from Minnesota Mining & Manufacturing Corp 
sit was further purified by evaporation at the tem- 
perature of a slush of dry ice in a 50-50 mixture of 
CHCI, and CCl, and subsequent condensation in a 
trap surrounded by liquid nitrogen, while pumping 
with a high-vacuum apparatus. Approximately 
385 g of final sample were obtained in this manner. 

\ mass-spectrographic analysis of the sample as 
used indicated a purity of 99.4 mole percent or 
better, with a maximum of 0.1 mole percent CO,, 
0.5 pereent No, and/or CO, and a trace of a hydro- 
carbon. A separate test with a phosphoric acid 
conductimetric method indicated only 0.02 mg of 
water per liter of sample. 


2.2. Apparatus and Method 


The construction of the flow calorimeter and its 
operation have been described in detail elsewhere 
(footnote 2); a brief résumé has also been published 
footnote 3). 

\ new electronic control system, designed and 
ult at the National Bureau of Standards, was used 


rt rk wes supported by the Ordnance Corps, Depertmert of the Army 
Masi and B. Petkof, J. Research N BS 48, 179 (1952). RP2303 
Masi, J. Am. Chem. Soc. 74, 4738 (1952) 


ace and J. G. Aston, J. Am. Chem. Soc. 70, 566 (1948 


successfully to control the temperature of the radia- 
tion shield to within +0.01 deg C of the average 
temperature of the tube containing the heated gas 
This control is to be described in a future publication 

The heat capacity was determined at a number of 
different flow rates at each temperature and pressure, 
so that any residual heat leak could be eliminated by 
extrapolation to the zero of reciprocal rate. “Blank” 
determinations were made at several rates at each 
temperature and pressure to measure the amount of 
cooling, 67, experienced by the gas in passing 
through the calorimeter when no heat was applied. 
A total of 67 heat-capacity experiments and 75 
blank experiments were made 


3. Results 
3.1. Heat Capacity of Hexafluoroethane 


If W is the power supplied to the heater, F is the 
rate of flow of gas, AT’ is the observed rise in temper- 
ature, and 67 is the fall in temperature in a corre- 
sponding blank experiment, the apparent heat ca- 
pacity is calculated by 

WFr-' 
AT +-6T 

The values of ( were corrected for the deviation 
of the observed mean pressure and mean temperature 
from the nominal values. These corrections were 
small in all The corrected heat capacity, 
called C, (observed), was plotted against the recipro- 
cal of the rate of flow at each temperature and pres- 
sure and straight lines were fitted to the points by 
the method of least squares. 

The least-square lines had very small slopes; more- 
over, the slopes were not constant with pressure at 
some of the temperatures. An analysis of variance 
was performed to test the significance of the slopes 
obtained. As a result, it was decided to abandon, 
in this case, the usual linear extrapolation to infinite 
rate, and to obtain the final values of heat capacity 
by simple averaging of the (, (observed) at each 
temperature and pressure. The results of this proc- 
ess are listed in table 1, converted to calories mole—' 
degree—'. (The calorie is 4.1840 j, the molecular 
weight is 138.02, and 0° C=273.16° K). Each of 
the experimental values listed in table 1 is the mean 
of at least four determinations. Following each tabu- 
lar value is the standard deviation of the experiments 
from the mean; the average of these deviations is 


(1) 


cases. 





TABLE | 


Tleat capacity of gaseous herafluoroethane 


Summary of results 


0.037 percent compared with + 0.036 percent for 
the corresponding precision index in the case of least- 
square lines fitted to the experiments 

The final values of heat capacity at finite pressures 
in table 1 were extrapolated linearly to zero pressure 
at each temperature, and the results are given as 
ay FAA observed.”’ It is believed that these heat- 
capacity values for the ideal gas C,F, are reliable to 
better than +0.15 percent 

An assignment of the fundamental frequencies of 
C,F, was made by Nielsen, Richards, and MeMurry.® 
Pace and Aston (footnote 4) have used this assign- 
ment in conjunction with their experimental third- 
law entropies at 176.61° and 194.87° K to obtain a 
value of 4,350 calories mole~! for the barrier to in- 
ternal rotation. A recent note by Mann and Plyler ° 
puts the lowest frequency (vy) at 220 cm™' by direct 
observation in the infrared, and suggests new esti- 
mates for the interatomic distances. The barrier en- 
ergy becomes 3,920 cal on the basis of these changes. 
The heat capacity has been calculated from this 
slightly revised assignment, for the harmonic-oscilla- 
tor, rigid-rotator approximation, and the results at 
the temperatures of the present experiments are given 
in table 1 as “C,°, spectroscopic.”’ They are seen to 
be 0.4 to 0.7 percent higher than the observed values 

3.2. Equation of State of Hexafluoroethane 

It is assumed here that in the low-pressure region 
covered by these experiments the equation of state 
of C.F, is of the form 

PV 

where B is a function of the 
thermodynamic relation 


( = ? ) — T( vs (3) 
oP /, oT" /p 
must be satisfied by the second derivative of the sec- 
ond virial coefficient B. The values of AC,/AP listed 
in the last line of table 1 were taken to be the left- 
hand member of eq (3), and B was assigned the form 
proposed by Hirschfelder, McClure, and Weeks: ? 
B=b—ce*’. (4) 
The values of ¢ and a were then chosen to fit the 
data of this experiment, whereas 6 was obtained 
from the vapor-density data of Pace and Aston 
(see footnote 4) at 25° C. The final expression for 


RT + BP, (2) 
The 


temperature. 


+ J. R. Nielsen, C. M. Richards, and H. L. McMurry, J. Chem. Phys. 16, 67 
(1048) i 
* 1D. E. Mann and E. K. Plyler. J. Chem. Phys. 21, 1116 (1953). 


'J. O. Hirsehfelder, E. T. McClure, and I. F 
201 (1942). 


Weeks, J. Chem. Phys. 10, 


10.00 
calories mol 


24. 718, +0. 011 
24. 617, +0. 08 
24. 526 
24.43 
24. 


+0. 010 


the virial coefficient in eq (2) was 


B91 


The experimental values of AC,/AP are plotted as 
as circles in figure 1. The calculated values obtaine 
by differentiating eq (5) are indicated by the soli 
curve, whereas the dashed curve is calculated fro 
the Berthelot equation, using the critical constants 
given by Pace and Aston (see footnote 4). Tabj 
2 shows that eq (5) reproduces the results of Pace an 
Aston’s calculation of B from vapor pressure and 
heat of vaporization measurements, at least as wy 
as the Berthelot equation. 


58.2¢/T em® mole 


, 


TaBLe 2 Comparison hetween equations of state 


he rafl uoroethane 


Value of 


Data of Pace 
ind Aston 
see foot 


note 4 


Berthelot 


Equatior 
equation scones 


cm mole 
S37 
SIS 

Ta 

76 

Al 

653 

611 

#24 

215 (ave 
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Effective Circuit Bandwidth for Noise with a 
Power-Law Spectrum 


Philip R. Karr 


The effective bandwith of tuned circuits for noise with a power-law spectrum is derived 


and discussed. 


noise bandwidth of tuned (RLC) circuits for 
white,” or “flat,’’ noise is well known.' However, 
orresponding results for noise with other types of 
ra have apparently not been discussed very 
ch. In connection with some noise-measurement 
work currently going on at the National Bureau of 
Standards, the writer worked out the effective noise 
bandwidth in the case of an f* spectrum (f being the 
frequency). The ratio between the effective band- 
width for the noise and the half-power band width is, in 
some cases, appreciably different from the well-known 
factor 1.57 that applies to flat noise; in other cases, 
the correction is quite small. These results are 
offered for publication in the belief that this knowl- 
edge will be of value to others engaged in noise work. 
Spectra of this form are encountered in transistor 
noise, resistor current noise, flicker effect noise, etc. 
We consider a random noise current whose spectral 
intensity (often called power spectrum) follows the 
law f*, driving a circuit consisting of the’ parallel ? 
elements R, L, C. 


The impedance Z of such a stage is 


Rd 


d+j( } - *) 


0 


vhere fo=(24LC)~'? is the resonance frequency, and 


Qnrf L l L 
RVC 


s the dissipation factor or reciprocal of Q. It is 
easily shown that df, is the frequency difference 
between the half-power points of the circuit; that 
is, it is the bandwidth 0. 

The mean-square output voltage of the circuit 
when driven by the noise is, with an obvious nor- 
malization, 


2rf,kRC YQ R 


- 


(V?) J, (Zi df (2) 


(2a) 


bd fek'l, 


Valley and Wallman, Vacuum tube amplifiers, p. 169 (McGraw-Hill 
", 1948) 


Inc., New York, N. Y 
its may be easily applied to the series RLC case 


Comparison is made with the usual case of ‘‘flat’’ noise. 


where 
reda 


dl? +( r 


The main task before us then is the evaluation 
of J. This evaluation is facilitated by the change 
of variable y=’, which gives 

l+a 
y 2 dy 


“+ 2y cos d 


where 


This integral may now easily be evaluated by con- 
tour integration in the complex plane. Omitting 
the details, the result. is 


if —3<a<l and —rZA< 4 

The restrictions on @ and X must be carefully 
adhered to. If valuesof a outside the specified range 
are used, the integral will diverge. Thus the device 
of assuming that the spectral law is constant for all 
frequencies, which appears to be justified in many 
cases when the letennal will converge, cannot be used 
when a>! or<—3. 


Now we may define bandwidth as 


the effective 
ie 
B= mp « 


so that 
B=bdI. (6) 
For flat noise we insert a=0 in (4), obtaining 
]=-2/2d, so that 


The 


which is the well-known result for flat noise. 
subscript refers to the case a=0. 





Using (7), (4), and (6), we may now write 


BB, is the correction factor for the effective notse 
bandwidth resulting from the deviation of a from 0 

The correction factor B/B, is an even function of 
a+1, that is for example, it is the same for a=0 as 
it is for a 2. (Therefore, we have the interesting 
result that the effective bandwidth for a 2 1s 
the same as it is for flat noise 

A plot of B/ By versus a for fixed d would take the 
form of a U-shaped curve with the bottom of the 
U occuring a 1: The curve is symmetrical 
about a B/By is less than 1 in the range 

2< a< 0), is equal to | for a 2 and a=0, and 
is greater than 1 outside this range; it increases 
toward infinity as @ approaches | from the left or 

3 from the right 

A family of curves for B/B, may be plotted without 
difficulty. However, a quick idea of the magnitude 


of the correction may be obtained from the fo! 
formulas to which (8) reduces in special cases 


B ) 
*). 


a 





For example, for d=0.1, a 1, we find immed 
ately B/B,~1—d/x =0.97. 

Thus for reasonably high values of Q(=1/d), «J 
corrections in the range —2<a<0 are rather smal 
Outside this range the corrections may be large 


Wasninecron, May 19, 1953. 
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Penetration of X- and Gamma Rays to Extremely 
Great Depths’ 


U. Fano ” 


Karlier work on the asymptotic trend of the X-ray intensity at great distances from a 
source is reviewed and completed in various aspects. The asymptotic law is shown to be the 
same as in the “‘straight-ahead” approximation (which disregards deflections) whether the 
primary energy is higher or lower than the energy of minimum absorption, provided a constant 
is replaced by the eigenvalue of a suitable Wick equation. The penetration in directions 
oblique to the source direction hardly ever attains its asymptotic trend when the source 
energy is lower than the energy of minimum absorption. This situation raises a difficult 
problem regarding t he penetration law in the range of great depths where the asvinptot 
trend is being approached very slowly 


l. Introduction 


The deep penetration of X-rays, or gamma rays, in an infinite, homogeneous medium 
are discussed. Much fundamental understanding of this phenomenon derives from the work 
of Wick [1]* on the analogous phenomenon of neutron penetration. Vice versa, the present 
work might contribute additional clarification to neutron problems as well as some background 
for still unsolved problems of charged-particle penetration. 

This paper constitutes a final report of developments since 1948. It includes a review of 
preliminary reports [2, 3, 4] and new material that was required to complete the initial program. 
Some concepts and techniques developed in the course of the study but not utilized in the even- 
tual solution are nevertheless reported briefly. 

The assumption of aa infinite, homogeneous medium disregards the effect of boundaries 
and inhomogeneities, which constitutes a separate, still largely unsolved, problem. 

The very deep penetration of X-rays depends primarily on the course of multiple Compton 
scattering under conditions where photoelectric absorption and pair production are compara- 
tively unimportant. At the low-energy end of the spectrum, for example, below 50 kev, where 
the energy shift of the scattered photons can be disregarded or treated as a small correction, 
the diffusion of photons has been studied by Chandrasekhar |5], but without specific reference 
to very deep penetrations. (Notice that photons below 50 kev disappear rapidly by photo- 
electric effect, in most materials.) At the high-energy end, where large amounts of X-rays are 
regenerated by electrons, a complete study of X-ray penetration would require a treatment of 
the whole cascade shower process. The X-ray regeneration remains moderate for secondary 
electrons up to 10 Mev in lead, aad up to 100 Mev in light materials, and may be treated, if 
necessary, in this energy range as a secondary source. The present paper disregards the re- 
generation of X-rays by electrons and treats pair production as a mechanism of outright absorp- 
tion. Nevertheless, its results apply to the penetration of the tail end of showers that is con- 
trolled by photons below 10 Mev in lead and below 100 Mev in light elements. 

We shall deal primarily with photons from 50 kev to 50 Mev which experience concurrent 
processes of energy degradation and multiple scattering terminating in outright absorption. 
Repeated Compton scattering of the “primary” radiation emitted by a source gives rise to 
“secondary’’ X-rays of lower energy traveling in all directions. In a medium of low atomic 
weight, a photon may be scattered 5 or 10 times, on the average, before eventual absorption 
by photoelectric effect. The main complication in our study arises from the generation of 


this complex radiation. 


Work supported by the Office of Naval Research and by the Atomic Energy Commission Reactor Division 
Appendix C is by L. V. Spencer 
Figures in brackets indicate the literature references at the end of this paper 
‘ This consideration applies also to the further development of Chandrasekhar’s method that was carried out by O’ Rourke [21] in connection 
with the X-ray penetration problem 
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This complication and the complicated dependence of the cross sections on the photon 
energy preclude a fully analytical treatment of the problem. Therefore, the analytical treat- 
ment that is developed in this paper does not lend itself readily to direct numerical applications. 
However, its results provide one of the footings for a practical method of numerical calcula- 
tion [6] whose range of application extends to very large distances from the source. The 
analytical laws of intensity variation versus depth also serve to extrapolate the results of num- 


erical calculations to still greater distances from the source. 

Throughout the present investigation, much reliance has been placed on gaining an 
initial qualitative understanding of the factors that control the process of degradation, scat- 
tering, and penetration of X-rays. Such an understanding enables one to adapt the mathe- 
matical procedure to the conditions affecting each specific problem. In this respect our work 
departs from other work in the same field that has been characterized by the initial adoption 
of less flexible procedures. For example, the X-ray distribution has been resolved into “orders 
of scattering” (that is, into once-, twice-, thrice-, . . ., scattered photons) [7, 8, 9]. The ex- 
pansion converges poorly where the photons experience protracted scattering and its applica- 
tion becomes cumbersome. However, at least in some instances, adequate calculations of 
penetration have been made [8]. Semiempirical improvements to the procedure have also 
been introduced [9]. Alternatively one may schematize the scattering and degradation 
process in such a way that the corresponding basic equations are sufficiently simplified to 
allow analytic solution {7, 9, 10, 11]. These simplifications, however, seem generally to in- 
volve considerable limitations in the applications of the theory and may delete some char- 
acteristic features of the phenomenon. Finally, calculations by the ‘““Montecarlo”’ sampling 
methods are possible [12] but do not yet appear convincingly successful, at least in their appli- 
cation to very deep penetrations, which requires the use of “biased sampling”’. 

Our treatment will be limited to X-ray distributions having a plane symmetry, that is, 
which vary in one direction of space only and are generated by sources having the same sym- 
metry. The study of spherically symmetrical distributions is essentially equivalent to that of 
plane-symmetrical ones. The relationship between these types of distributions is well known 
[1] and will be indicated in Appendix A. We shall also assume that the source is concentrated 
on a single plane, which implies no additional restriction owing to the linearity of the problem. 
Notice that plane symmetry may always be attained by considering, instead of the X-ray 
flux at each point and in each direction, only the integral or the average flux over all points 
of a plane and over all directions forming the same angle with the plane. The reason is that 
this plane integral flux generated by a localized source is equal to the localized flux generated 
bv a plane distribution of sources. 

The most general type of source may be regarded as an aggregate of point monodirectional 
sources. The study of X-ray penetration under conditions of plane symmetry may serve as a 
first step for the treatment of point monodirectional sources. This approach has been ex- 
ploited in the study of X-ray penetration to small or moderate distances from a source [13]. 


2. Survey of the Problem 


This section contains a qualitative analysis of the deep penetration of X-rays, with the 
purpose of pointing out the specific problems that require a detailed mathematical treatment. 


2.1. Transient Processes and Equilibrium States 


The layers of material near an X-ray source are traversed almost exclusively by primary 
radiation because it takes a certain thickness of material to generate a substantial amount of 
secondary radiation. Therefore, one expects the intensity of secondary radiation to build 
up rapidly from one layer to the next near the source in a sort of transient process. Vice versa, 
at greater depth within a material one may hope to find some sort of steady state in which the 
secondary radiation is present in substantial amounts. The nature of this steady state is the 
main object of this investigation. 





A typical steady state in which a softer radiation continually arises from a harder one is 
ommonly called a state of radiative “equilibrium.”’” The X-ray penetration problem does not 
ead to a typical state of equilibrium, but it is helpful to review first the properties of such a 
state. Equilibrium involves, in general, three different features, namely: (1) The ratio of the 


“ce 


ntensity of the secondary to the intensity of the primary radiation approaches asymptoticall 
a maximum value, as the depth of penetration increases. (2) The quality of the secondary 
radiation also becomes independent of the depth of penetration, because the intensities of its 
various spectral components bear a constant asymptotic ratio to the primary intensity. (3) 
Therefore, the intensity of the whole radiation is controlled by the progressive attenuation of 
the primary radiation, which follows an exponential law in the case of X-rays. If these 
circumstances applied, they would greatly simplify the study of penetration. The maximum 
values of intensity ratios mentioned in (1) and (2) are determined by the ratios between the 
attenuation coefficients ° « of the X-rays of various energies. The depth of penetration at 
which the intensity ratio between any two X-ray components attains a value within, say 10 


percent of the maximum, turns out to be inversely related to the difference between the attenu- 


ation coefficients of those components. 


2.2. State of Limited Equilibrium 


The approach to a steady state under the conditions of hard X-ray penetration presents 
the following characteristic departure from the more familiar equilibrium conditions. In the 
Compton effect, some of the X-rays have very nearly the same energy, the same direction, and 
the same attenuation coefficient after as before the scattering. In fact, the attenuation coeffi- 
cient of some among the secondary X-ray components differs only infinitesimally from the 
attenuation coefficient of the most penetrating X-rays that are present. This circumstance 
does not merely slow down the approach to equilibrium, it actually suppresses the features (1 
and (3) in that the intensity ratio of secondary X-ray components to a monochromatic primary 
may grow beyond any limit. 

At the same time one may expect that, within a finite thickness of shield, some sort of 
equilibrium should be attained among those components of the secondary X-ray spectrum, 
whose attenuation coefficient differs from that of the most penetrating component by a suffi- 
ciently large amount. In fact, the relative intensity of most secondary components and, 
therefore, the average ‘‘quality”’ of most of the secondary radiation turn out to approach a 
limiting value within the thickness of a finite shield, as indicated in feature (2). On the other 
hand, features (1) and (3) do not obtain. The intensity ratio of this secondary radiation to 
the primary X-rays does not approach a maximum value, but it keeps increasing; the attenua- 
tion of the primary radiation alone does not control the attenuation of the secondary radiation. 

In a hypothetical state of true equilibrium, the flow of radiation at great depths within a 
shield is described as the product of an exponential function of the depth and of a function of 
the energy and direction of each secondary radiation component. In the sort of steady state 
that is actually achieved, the flow of most of the secondary radiation at great depths is described 
as a product of a nonerponential function of the depth and of a function of the energy and 
direction of each component. Calling z the depth of penetration, measured from the source, 
FE the energy of a photon, and @ the angle that its direction of travel forms with the + axis, 
ky the energy of the primary photons, the flux of photons of various energies and directions 
(7,40) takes then the form 

d(x, FE, 0)~ f(x) g(E, 0) (1) 


when z is sufficiently large and £ sufficiently smaller than £5. 

At large depths most of the energy is carried by the softer secondary components. One 
ean, therefore, use the approximate form of the flux to evaluate the quantities of practical 
interest. These quantities vary, as a function of depth, in proportion to f(z). 


The term “attenuation coefficient” is used instead of “‘absorption coefficient” to make it clearly understood that » includes the effect of 
moval of X-rays from a narrow beam by Compton effect in addition to their true absorption by photoelectric effect or pair production 
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2.3. Intensity Variation at Great Depth 


The over-all attenuation of the radiation does not follow an exponential law because the 


most penetrating secondary components never approach an equilibrium. The determination 


of the over-all course of attenuation, which is described by the function f(x), requires, therefore, 


a study of the ever developing transient process of formation and destruction of the hardest 


secondaries. This study constitutes the main topic of the present paper. 

By centering one’s attention on the formation and absorption of the most penetrating 
secondaries, one attains a substantial simplification of the problem. The secondaries that 
need be considered cover only a narrow spectral range. Therefore, the variations of the 
and absorption cross sections within this spectral range can be treated as small 


scattering 
Furthermore, the greatest contribution to the deep penetration arises from 


quantities 
secondaries directed in a very narrow beam perpendicular to the source plane.® 

However, the small differences of attenuation coefficient and of direction among the most 
penetrating components cannot be disregarded, even in the first approximation. On the 
contrary, it is just these differences that control the relative intensity of the various components 
and thereby the over-all course of attenuation at great depths. The following argument indicates 
that small changes of attenuation coefficient and of direction have comparable importance 
Radiation of wavelength A, which travels in a direction forming an angle 38 with the axis of 
penetration, s, travels a distance J/cos J as its distance from the source plane increases from 
rtor+aA. Therefore, its intensity decays along zs with an effective attenuation coefficient 
Compton seattering of a photon, with a wavelength shift of 6s Compton units, 
results in a change of both w(A) and cos 8. Ifa photon starts with a wavelength \ in a direction 
|, after one scattering its directioa changes, according to the Compton 
law, to 8—arecos 1] 5A] and its effective absorption coefficient from (Ay) /Cosdy= u(y) to 
u(Ay + BA) [LE — SA] ~ wl Ay) + (du dd) SA+ w(AQ)SA. The corrective term (du/dd)dX arises from the 
change of M, the term u(A)bA from the change of cos 8. The ratio of the two corrective terms 
is d log w/dd, with X in Compton units, which is a number of the order of 1 for most materials 


ul(A) cos J 


with @—0, cos 6, 


and wavelengths of interest 
The importance of the small deflections was not appreciated in the early stages of the 


investigation, which disregarded altogether the changes of direction experienced by the most 
penetrating radiation components, that is, relied on the “straight-ahead approximation.” 
The results obtained in this manner, which are reported in [2, 3] and reviewed in section 5, 
proved useful, nevertheless, because the effect of small deflections does not change the analytical 
form of the intensity variation but only the value of certain numerical constants. This simple 
mode of action of the deflections emerges as a result of the analytical treatment of section 7 
and may be visualized as follows. Photons that travel in directions oblique to the direction 
of maximum penetration are selected-against in the course of deep penetration, much as if 
they had, in effect, a larger attenuation coefficient (A) (see above). The selection operates 
simultaneously in favor of photons with near-minimum yu and near-minimum, that is, near-zero, 
obliquity 8 As a result, photons of each wavelength are concentrated in a cone of directions 
whose aperture is the narrower the nearer to the minimum is «. This is to say that the effect 
of small obliquity merely parallels and amplifies the effect of small variations of the attenuation 
coefficient. Indeed, the two effects result in a fixed ratio to one another (see sections 6 and 
7), so that the study of the effect of variations of » alone, in the straight-ahead approximation, 
gives the correct analytical form for the trend of the total intensity at great depths. 

Two different analytical forms for this trend result under different circumstances, depending 
on the type of variation of the attenuation coefficient among the most penetrating components 
As long as the energy of the primary X-rays is not too large, none of the secondary X-rays 
is more penetrating than the primary X-rays. The hardest secondaries are then simply those 
whose energy is just a little lower than the energy of the primaries. The progressive accumula- 


* When the whole primary radiation travels in directions oblique to the source plane, there occurs a more complicated situation that will be 


discussed in section 2.6 and in section 8 





on of secondary components under these conditions has been studied in |2] under the 
raight-ahead approximation and in [4], taking into account the deflection effects. The 
tensity distribution law f(r) is shown there, and in section 9, to be of the form 


fir) rx exp( pot). (2) 
lere wo is the attenuation coefficient of the primary X-rays and A is a number, usually of the 
wder of 1; K depends on the cross section for Compton scattering and on the first derivative 
f the attenuation coefficient with respect to energy, evaluated at the energy of the primary 
X-ravs. (For numerical values of A and the parameters on which it depends, see appendix 
C, tables 1 and 2.) The exponential factor represents the usual effect of attenuation of the 
primaries, Whereas the factor «* represents the effect of accumulation of hard secondaries and 
is called the “buildup factor’. Numerical values of A are given in appendix C 

At very high energies, absorption by pair production mecreases So rapidly that higher 
energy X-rays are actually less penetrating than lower energy ones. Therefore, in the case of 
high-energy primaries, some of the secondary X-rays are more penetrating than the primaries 
The attenuation of the total radiation depends upon the formation and decay of the most 
penetrating secondary components, whose photon energies may be much lower than that of 
the primaries (in lead the hardest photons have about 3-Mev energy). This phenomenon has 
been treated im [3] in the straight-ahead approximation, but the effect of small deflections 
has not been taken into account in previous papers. Both treatments are given in sections 5, 
7, and 9 and show the intensity distribution law f(z) to be of the form 


, exp |/7 (u» 2) jexp mt). 


Here yu, Indicates the minimum value of the attenuation coefficient, that is, the attenuation 
coefficient of the most penetrating X-rays in the material under consideration. The factor 
exp M»t) represents the attenuation of these X-rays The factor 27° exp|/7(u»,2"] is the 
buildup factor, corresponding to 2* in (2), and // is a constant that depends on the cross section 
for Compton scattering and on the values of the attenuation coefficient and of its second 
derivative, evaluated at the energy of the most penetrating X-rays. (For values of u,, and I, 


see appendix C, table 3 
2.4. Numerical Methods Combined With Analytical Results 


The analytical results, (2) and (3), represent limiting forms of the intensity distribution, 


valid at extremely great depths, where the distribution depends on the generation and attenua- 
tion of an extremely narrow range of secondary spectral components. Those results must be 


supplemented in two directions. In the first place, one must find not only the relative X-ray 


intensity at different positions but the absolute intensity, particularly of lower energy compo- 
nents traveling in various directions, that is, the function g(/,3) of eq (1). In the second place, 
one must be able to carry out calculations without excessively stringent limitations on the 
important range of secondary spectral components. 

Various partially successful attempts were made to develop methods of successive approxi- 
mation to take into account successively higher derivatives of the cross sections in the important 
spectral range. However, much greater success was eventually achieved by the semiasvmptotic 
numerical method of Spencer [6], which fulfills both of tne requirements indicated above 
Additional numerical results for the X-ray distribution at moderate depths are provided by 
polynomial method of calculation [13]. 

Reliance on numerical methods to the extent of solving the transport equation numerically 
was dictated by the complicated dependence of the cross sections upon the X-ray energy. At 
the same time, since the transport equation involves three independent variables (namely 
distances from the source, photon energy, and direction of propagation), the reduction of the 
iumerical burden to manageable proportion requires much guidance from qualitative analysis 
ind as much help from analytical development as conveniently possible. Spencer has stressed 





the rapid gain in computational efficiency that can be derived from general information on th, 


desired solution [6, 14]. Economy also arises from suitable substitutions of the variables of 


the problem (see, for example [13]). 

The space variable, depth of penetration from the source, can be separated from the othe: 
variables by a Fourier-Laplace transform. This transformation has the additional advantag: 
of replacing the depth of penetration, which varies over a wide range, with a transform variable 
whose important range is narrow. As will be seen in section 4, the analytical treatment that 
determines the intensity distribution laws (2) and (3) deals specifically with the behavior 
of the distribution function at a singular point of the transform variable. The numerical! 
method holds only for values of the transform variable at a finite distance from the singular 
point. The analytical treatment serves as an essential complement to determine the critical 
behavior at infinitesimal distances from the singularity. Furthermore, the knowledge of this 
analytical behavior and of the resulting behavior of the corresponding inverse transform, 
(the function f(z) of eq (1), (2), and (3)), serves as a guide for the inverse transform procedure 
that is also required in the numerical work. 

The two remaining variables, photon energy and direction of propagation, can be separated 
oaly within the limits of application of the analytical treatment, that is, for extremely deep 
penetration and narrow range of photon energies (section 7). Otherwise, the interlinkage of 
energy and direction offers the most serious difficulty to the solution of the problem. The 
difficulty has been overcome in the numerical work by learning how to take into account this 
interlinkage without excessive complication [6]. This is done by describing the directional 
distribution adequately by means of a few parameters; namely, moments when the distribution 
is peaked, and coefficients of a Legendre expansion when it is flat. The interlinkage of energy 
and direction is reviewed briefly in section 6. Numerical methods have to complement the 
analytical treatment even when the separation of variables succeeds in order to work out the 


directional distribution (section 7 and appendix B). 
2.5. Verification of the Qualitative Analysis 


The introductory qualitative picture of the X-ray penetration, which has led to the results 
indicated by eq (1), (2) and (3), has been verified to a considerable extent by experimentation 
and by independent numerical calculations. 

Two experiments dealt with a point source of Co™y-rays surrounded by a mass of water 
In the first experiment [15] ionization-chamber and Geiger-counter measurements of total 
intensity were made at various distances from the source. The results of these measurements 


should be represented, in terms of the approximate eq (1), by f(z) 27 sin ode | g(Ed)r(E\dk, 


J0 0 
where r(/) indicates the response of the measuring instrument. In the absence of information 
on g(/,3), the plot of the measurement versus the depth of penetration z should follow the trend of 
the function f(z). The plot does, in fact, take up the trend of (2) at great distance from the 
source (1.5 to 2.5 m of water), with a value of K equal to that predicted by the detailed theory of 
[4] and of section 7. 

In the second experiment [16], the spectrum of the secondary electrons in the water was 
measured at various distances from the source. The plot of the number of electrons of various 
"Oe sin odd [ g(E,8)R (E,edE, 


energies ¢, should be represented according to eq (1) by f(z) 
0 J0 


where 2(E,e) indicates the probability of production of an electron of energy ¢ by a photon of 
energy &. Qualitatively, the shape of the plot becomes independent of the penetration depth 
z at great depths as predicted by (1). 

Equation (1) gives no information about the shape or the absolute magnitude of g(E,%). 
However, rather detailed calculations of the complete distribution function #(z,£,8), for all 
depths of penetration covered by the experiments can be made by the independent method of 
polynomial expansion [13]. These calculations have successfully predicted the quantitative 
results of the two experiments. 
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In addition, the results of these calculations for moderately great depths provide a semi- 
quantitative check of the predictions of the theory of penetration to very great depths. These 
predictions find good verification. In particular, the trend toward an equilibrium distribution 
for the lower-energy secondary X-rays is clearly displayed both with respect to the spectral [13] 
and to the directional 17] distributions (see fig. 1). The emergence of an equilibrium state at 
lower energies, while the higher energy components persist in a state of nonequilibrium is 


most clearly displayed in the results of semiasymptotic calculations [6] (see fig. 2 


2.6. Influence of the Source Direction on the Penetration Law 


Previous theoretical work on the penetration of radiation from plane sources has dealt 


only with isotropic sources and with sources concentrated in a direction perpendicular to the 


source plane. In general, one may want to deal with source radiation aimed in directions 
that form a fixed skew angle with the source plane. However, this new problem raises questions 
that are largely unsolved and that will be indicated briefly here and discussed somewhat further 
In section S 

At the outset one must distinguish between the conditions where the energy of the primary 
X-rays lies below and above the energy of minimum attenuation (corresponding to eq (2) and 
(3), respectively In the first event the penetration law is controlled by secondary X-rays of 
energies just below the source energy. The overwhelming majority of these secondary X-rays 
travel in directions very nearly as oblique as those of the primary X-rays. Therefore, the 
obliquity of the source has a very major influence on the intensity distribution at very great 
depths 

In the second case, when the source energy exceeds the energy of the most penetrating 
secondaries by a substantial amount, the obliquity of these secondaries is less immediately 
related to the obliquity of the primaries. It will be shown in section 8 that the penetration is 
controlled by the accumulation of secondaries that have been degraded by steps to the energy 
of minimum attenuation and simultaneously deflected to directions nearly perpendicular to 
the source plane. The absolute intensity of these secondaries depends, of course, on the 
obliquity of the source but their rate of accumulation and, hence, the law of deep penetration 
are independent thereof 

Therefore, the main difficulties arising from source obliquity occur only when the energy 
of the primaries lies below the energy of minimum attenuation. If % indicates the angle 
formed by the source direction with the z direction, perpendicular to the source plane, the 
“effective attenuation coefficient” of the primary X-rays of wave length Ay is w(Ao) /cos dB) (see 
section 2.3 After one Compton scattering with a wavelength shift 6A, the obliquity @ is 
given by cos J =cos J)| | — 4A] +sin J[25A—4d*]'? cos g, Where g may have any value. Minimum 
effective attenuation after one scattering is obtained when g=0, but this minimum value 
u(Ao + 6A)/< cos J >a, may be larger or smaller than the initial value u(y) /cos dj), depending on 
the values of du/d\ and of d&). The most favorable reduction of the effective attenuation co- 
efficient is brought about by a succession of scatterings with infinitesimal 4\’s, all directed so 
as to reduce 3. An infinite sequence of such scatterings would eventually bring about a deflee- 
tion all the way to 3-0 without any finite increase of u, but the probability of such an event 
is, of course, infinitesimal of a high order. Nevertheless, an unlikely process of this kind must 
determine the penetration at extremely great depths. That is the penetration will be eventually 
controlled by the seanty accumulation of secondary X-rays with energies very near the pri- 
maries and direction nearly perpendicular to the source plane, that is, with effective attenuation 
coefficients near wo. The course of accumulation of secondary X-rays that are effectively more 
penetrating than the primaries because of reduced obliquity has not yet been studied. | Knowl- 
edge of the limiting trend at extremely great depths may not be very important because it may 
not even become recognizable except under unrealistic conditions. Nevertheless, lack of 
understanding of this process may undermine the effectiveness of the semiasymptotic and 
polynomial methods of calculation, whose application is now being planned. 
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If the source emits X-rays of various obliquities, some of which travel initially in the 
rection perpendicular to the source plane, it will be just the X-rays in this direction and 































eir secondaries that achieve great penetration. Deeper and deeper penetrations arise from 
rrower and narrower bundles of primary X-rays. Therefore, a source with a broad distribu- 
yn of directions contributes less and less toward deep penetration, as compared to a source 
med exactly in the direction perpendicular to the source (d)=0). This argument will be 
mfirmed by the analytical treatment of section 8. In particular an isotropic source, or any 
ource with practically uniform distribution in directions nearly perpenducilar to the source 
Jane, vields an intensity distribution law in accordance with eq (2) but with a value of AK 
ower by 1 than the value for a source concentrated at d)=0 [6] 
3. Transport Equation 
The following symbols will be used: 
r=(x,y,z)=coordinates of a point of a medium 
A\=wavelength of a photon in Compton units (h/m« 
u Ur Uy Ue v.¢o unit vector indicating a photon direction 
u(r total (narrow-beam) attenuation coefficient of the medium for photons of wave- 
length 
ur» = Thomson scattering coefficient (probability per unit path) of the medium 
low energy limit of the integral Klein-Nishina coefficient 


, 


aN 3/8) ural dX’ /A+A/d N’/d)7(A— A’ — 2) (A—X’)] = Klein-Nishina differential 
scattering coefficient of the medium for Compton scattering with a wavelength 
change from \’ to A, per unit A, for A’ <A<A’ +2. 

Y(r.u,4) number of photons per unit volume, per unit solid angle and per unit 


Og Dirac’s delta function 

\ wavelength of the source photons (if monochromatic 

Lb Mh AN , 

A, = wavelength of the photons whose attenuation coefficient has the lowest value 

} in the range of integration from \, to X 
\, =Wwavelength of the photons whose attenuation coefficient is the absolute mini- 
mum in the medium under consideration. 
ben = pA, 


The degradation, penetration, and diffusion of X-rays is governed by a transport equation 
his equation represents the rate of change of the density of photons with a specified direction 
ind wavelength, from one point to the next in the direction of propagation, as the sum of three 
terms namely: (a) the attenuation of the photon density as a result of absorption and scattering, 
b) the addition of photons that take up the specified direction and wavelength as a result of 
Compton scattering, (c) the addition of photons with the specified direction and wavelength 


from the source. The equation is: 


u-grad ) (r,u,\ u(A) V(r,u,r 


on s - 
dn’k(n’,d du’ (2r)'6(1—u-u’ —A+ ND’) Yir,u’,r S(r,u,), (4) 
2 ir 


. - 


vhere S(r,u, \) represents the source, that is, the density of photons produced at r with the 

lirection u and the wavelength \. The 6-function represents the Compton law which requires 
to equal \’+ (l—wu’ 

Following the plan indicated in the introduction, we assume that the source distribution 

is a plane symmetry, that is, that it does not depend on y and z, nor on the component ¢ of u, 


uit only on d=arceos u,). Asa result, Y will be similarly independent of y,z and ¢g. We also 
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assume that the source is concentrated on the plane s=0 and at a wavelength }\=o, but this 
constitutes no real restriction owing to the linearity of the equation. Accordingly, we have 


S(r.a.d) = 6(7)5(A— Xo) f(u,) and the equation reduces to 


u,OY (r,u,,r)/ Oz u(A) VY (x,u,,d) 
AN kiN Xd du’ (2r)~'6(1—w-u’ — A+ WN’) V (4,0), ’) +-6(2)6(A— &) f(u,). (5) 
2 Jin 


Notice that the kernel &(\’,A) of this equation may be multiplied by any ratio g(\’)/g(\), 


provided Y is suitably renormalized 


4. Fourier-Laplace Transform 


The transform method involves a separation of variables. It represents the X-ray distri- 
bution as a superposition of components, each of which has a certain angular and spectral dis- 
tribution uniform over all the space. The intensity of each component varies exponentially or 
sinusoidally from point to point according to exp(—pzr), where p may be complex.’ 

The transform and its inverse are: 


y(pyuz,r exp (pr Y(r,u,,r)dz, 


- 


, | re exp pr y PyUzd dp. 
J—i« 


® The transform method of analysis consists of two steps: (a) A study of the directional and 
spectral distribution y(p,u,,d) for various values of p, (b) an evaluation of the inverse transform 
integral (7). 

The evaluation of the integral may proceed by the path of steepest descent. Whether this 
method is eventually followed or not, an analysis of the “‘topography”’ of the transform y(p,u,,d), 
that is, of the distribution of its absolute values over the complex plane p proves very useful. 
This analysis indicates what region of the plane yields the largest contribution to the integral 
when the steepest descent method is followed. The determination of y(p,u,,d) in this region, 
which is often very limited, proves anyhow to be very efficient for evaluating the inverse 
transform. 

The equation that governs the directional and spectral distribution y(p,u,, ), for each 
value of p, derives from (5) by means of the transformation (6), that is, by multiplication by 


exp(pr) followed by integration. This equation is: 


[u(A) — pu,ly(pyuz,r) | dw kX») | du’ (2x)~'3(1—u-u— d+ d’)y(p,u,, dX’) +5(A— &) f(u,). (8) 
0 4x 

The source term in this equation does not depend on p because the source term in (5) is con- 

centrated on the plane z=0. 

Because the wavelength of a photon increases as a result of Compton scattering, the dis- 
tribution of photons of each wavelength \ depends on the distribution of photons of shorter 
wavelengths \’ but not on that of longer wavelength photons. If one knows the photon dis- 
tribution (in spectrum and direction) at all wavelengths \’<\, one can evaluate the integral 
in (8) and thereby find the distribution y(p,u,, 4) on the left side of the equation. Thus eq (8) 
may be solved in principle stepwise, proceeding from the source wavelength \) to longer and 
longer wavelengths. This procedure, which reflects the course of the physical process of energy 
degradation, points to an important effect. If the intensity of X-rays of a certain wavelength 
builds up to a high value, Compton scattering of these X-rays builds up, in turn, the intensity 
of all longer wavelength radiation. 


' Notice that p has here a sign opposite to that which is most common in the literature. The convention adopted here proves convenient in 


problems of penetration and straggling 
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rE 3. Elements of a relief map diagram of the trans- 


. “IGURE ) r ¢ , / th l 
, distribution yon the plane oj the complex variable Pp. Fe a 1. y tata oJ ; plot of iy P Coney oS VERS Cues 


of figure 
Lines over which y(p) is singular trend of the steepest 
nt line of y(p) exp (—pr) for positive p , approximate position 
saddle-point 
Equation (8) shows that y(p,u,, \) becomes infinitely large when u(A)— pu, vanishes. Fur- 
thermore, ¥ grows rapidly for those combinations of p,v, and \ for which u(A)— pu, is particu- 
larly small. These combinations of variables are critical in various respects. We just men- 
tioned that once y(p,u,, \) has become very large for certain values of uw, and \, it remains 
correspondingly large for all longer wavelengths. We also know that the distribution of large 
values of y(p,u,, ) determines the path of steepest descent of the integral (7). Finally, the 
process of numerical integration becomes increasingly difficult as the rate of rise of y(p,u,, 
increases, so that one may want to utilize, at least locally, analytical methods of solution. 
The coefficient u(A)— pu, of (8) vanishes only for positive real values of Pp larger than u(\) 
or for negative real values smaller than —y(d). Therefore, eq (8) becomes singular only for 
values of the complex variable p which are confined to the sections of the real axis marked 
heavily in figure 3. The lines terminate at p=-+4y,, where yu, indicates the smallest value of 
u(X’) in the range of wavelengths from \) to’. There appears to be no other point in the plane 
p where y(p,u;, d) diverges 
Figure 4 shows a schematic plot of |y against real values of p. The plot rises to infinity 
at p= +u,, where the equation has become singular for u,=1 and for the wavelength A, <) at 
which uw takes its smallest value u, 
The integrand of (7) equals y(p,u,, \) times exp(—pr). The factor exp(— pr), regarded as 
a function of p, slopes down to the right for positive values of x and to the left for negative 
values. The slope increases in proportion to the distance from the source, 7. Therefore 
the integrand has a minimum, for real p, where 0 log y(p,u,A)/Op=~s, that is, at a point that 
approaches p=, when z is large and positive, or p u, When sz is large and negative 
This lowest point of exp(— pr)y(p,u,d) against real values of p constitutes the one and only 
saddle point along the steepest descent path when exp(— pzr)y(p,v,, ) is integrated with respect 
) to the complex variable p (see fig. 3 
The minimum and the saddle are sharp when |z| is large. Under these conditions, a knowl- 
edge of y(p,u,, \) over a narrow range of real values of p suffices to evaluate the integral (7) 
and thus to calculate the total distribution of X-rays far from the source. The main object 
of this paper is just to study the trend of y for real values of p near wu, (or —y,). Section 9 
deals with the actual evaluation of the integral (7). 
For smaller values of the depth of penetration z, the saddle point lies nearer to p=0. 
The direct calculation of Y(x,u,,4) at moderate depths also provides convenient expansions 
of y in powers of p, or of related variables like p/(u,—p) and p?/(u;—p*) [13]. This infor- 
mation may serve to interpolate y between the proximity of p=0 and of p= + uy. 
5. Special Analytical Solutions—Straight-ahead Approximation 
In order to study the trend of y for real values of p near u,, we seek approximate analytical 
olutions of the transform equation (8) valid for values of \X~\, and of u,~1. This search 


vill be conducted first in the straight-ahead approximation even though this approximation is 
adequate. 
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The straight-ahead approximation consists of setting u,=1 on the left side of (5) or (8 
This procedure is equivalent to the replacement of the depth of penetration z in (5) with the 
path length traveled by the photons, irrespective of the obliquity of penetration [18]. The 
variables uw, and uw that still remain in the eq (8) can then be eliminated by integrating over u 


setting 


| duy(p,uz,d) ( pyr 
Jae 
we obtain the straight-ahead approximation equation 


([u(X) — p]®(p,d) k(’,d)(p,d’)dd! + 8K — Ag), 


where I f(u,)du 


oie 
Because we are interested in an approximate solution valid over a narrow range of variation 
of \, we disregard the variations of the differential scattering coefficient & over an interval 
\’~d and take 
kN) RO A) = (3/4) umm =C, (11) 
in agreement with the definition of & (see section 3). 

Before entering this approximation into (10), we must distinguish two different situations, as 
in section 2.3. (a) If the primary radiation is the most penetrating one, that is, if ¢,= (Ao) = wo, 
the range of approximate integration over \ is right next to the source wavelength ). 
(b) If the most penetrating radiation is that of wavelength \,,, for which uv has its absolute mini- 
mum value u,,=u(d,), that is, if <A, the range of approximate integration cannot run all 
the way from \y to A. In this event, we do not attempt to solve the inhomogeneous equa- 
tion (10), but only the corresponding homogeneous equation, in the range of \~),.. 

Therefore, we seek the solution of 


A 
[u(X PIPO,p) Cfo, pan’ -T 6(A—r% 


A~ Ny > Am, 


7» 
and of pI®O,p) c| @(d’,pidn’ 
for 
A~ Am > Do- (12b) 


These equations reduce to differential form, if a derivative is taken with respect to A, and 
their solutions are 


mr 
P(A\,p) = 1T(wy— p)'C [w(A) — p]“ exp iC [_ an’ [u(r’) pit (13a 


7A 
¥(A,p) =C [u(A)—p)"' exp CJ dn’ /[w(X’) pit (13b) 

These solutions are derived as special cases of a more general procedure in {18 section 11b] 

We are interested in the behavior of these solutions when p approaches the singular point 
uy, that is, w in (13a) and g,, in (13b). 

In case (a) the first derivative of u(d’) at the singular point yo, (du/dd’)o= go, is different 
from zero. Therefore, the integral in the exponent of (13a) has a logarithmic singularity. 
In case (b) the first derivative of u(A\’) vanishes at u(A’)=u,, ’=A~. Therefore the integral 
in the exponent of (13b) over an interval of \’ ~,, is of the arctangent type and its singularity 
has the form (u,—p)~'”. If u(d’) is regarded as linear in case (a) and parabolic in case (b), 
the solutions (12) take the characteristic singular forms [3]: 





C 1() 
P| exp }= log © Cin “o-' (14a) 
(uo Ho 


P(A,p) ~C [w(dA) — pp]! exp { 22C/{: rt. (14b) 


where pb» = (d7y/dX*)xm.- 

Even though the solution of the transform equation (10) for values of \ A, cannot be 
»btained in this simple manner, its dependence on p, as p approaches y,, is described by (14) 
for the following reason. As soon as X is sufficiently larger than \, so that (A) — p Me—P, 
p may be replaced in the equation by u,,.. This implies that the further variations of ©(\,p), 
as \ keeps increasing, are independent of p, that is, that a state of limited equilibrium, within 
the meaning of section 2.2, has been achieved. As p approaches u,, \ need depart from \, less 
and less, before the condition u(A)— p us — pis adequately fulfilled. Therefore the approxi- 
mations made in the derivation become increasingly good. 

The inversion of Laplace transforms with the singularities (14a) and (14b) will be con- 
sidered in section 9. However, we anticipate here that the inverse transforms behave, respec- 
tively, as (2) and (3) for large depths of penetration z. 

For the sake of completeness, we mention yet another situation, which combines the 
features of (a) and (b), namely, the case where \y<A<A,,,A,=A. Here the pertinent form of 
the equation is (12b) and its solution (13b). However, the derivative (du/dd), at A=, is 
different from zero. Hence, the integral in the exponent of (13b) has a logarithmic singularity 
and we find, instead of (14b), 

O(A,p) ~€ [w(A) — pO ee (14e) 


Various attempts were made in the past to develop methods of successive approximations, 
starting from the approximate solutions (13) and (14), to solve the straight-ahead equation (9) 
for all values of A. Direct numerical integration eventually proved superior [6], even for values 
of p rather near y,, despite the fact that the rapid rise of y along \, for \ near \,, requires small 
steps of numerical integration. 


6. Possible Treatments of the Obliquity of Penetration 


As pointed out in section 2.3, the deep penetration of X-rays is influenced to a substantial 


extent by the small obliquity of the most penetrating components. In order to evaluate this 
effect, one must, in principle, work out the directional distribution of the X-rays of each wave- 
length, that is, solve the complete two-variable equation (8). In practice, one may want to 
learn about the distribution-in-angle only the least amount necessary to evaluate the effect on 
the trend of penetration 

In order to analyze this situation, one may expand the distribution-in-angle of y(p,u,,d) 
into Legendre polynominals P,(u,), 


y(pyuz,r) > (2 L-1)(49r) ‘yi (p,®) Pi(u,) (15) 
with 


Yi(p,X) J duP (u,)y(pjuz,). (16) 
: dn . 
Equation (6) multiplied by each of the P,’s and integrated, reduces to the system of equations 


uA) Yo(p, >») — py (P,®) [PFO mr )an +60 do) J, f(us)d 


u(A)ye(p,r) — pL + Dy yiga (pd) + lyri (p,d))/ (21+ 1) 


‘. k(N A) PL — A+ WM’) yp dN +8(.A— & i” P(u.f(u,jdu. 
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The first equation of this system has the same form as the straight-ahead equation (10), 
except for the fact that y, appears in the second term on the left side instead of yo. If there 


were a way of guessing the ratio 
Yi (Pr) /Yo(Py® f, u,y(p,uy,r)du ( y(p,u,,r)du u,(p,r), (18) 
: ‘ “hie Jax 


y, could be replaced with 7,(p,A)yo(p,d) in the first eq. (17), which could then be solved as a 
single-variable-straight-ahead equation. Thus a knowledge of @,(p,d), or of an equivalent 
parameter, embodies all the information about the directional distribution that is immediately 
relevant to the penetration problem. 

Instead of u, one may want to use the equivalent parameter [6] 


gi(p,¥) =1—U,(p,d) (l—u,)y(p,u,,r)du y(p,u,)du, (19) 


Jinx 


which represents the first moment of the angular distribution. If we set y,;=(1 91) Yo in (17), 


the first equation becomes 


*) . 
[u(r p- pa pyr ly (p,r kiN’ Yo( pv dn +-6(X X | f u, du. (20) 
J 0 de 


Thus, g;(p,d) is seen to measure the departure of the eq (20) from the straight-ahead eq (10). 
In the critical range of wavelengths, \~.,,the solution of(20) may be expressed in a manner 
similar, for example, to (13b), that is, 
we aT | 
Yo(p,yr ~f [u(r Pp pa p,r)] exp) dx [u(r ) p- Pai( pyr IK 
Cl w(r) p| : (> Cc dy’ ) 


, (21) 
Pair, p)/[w(r) p| “*P) 1+ pai(p.r’)/[w)—p] wr) —ps 


The last expression differs from (13b) in that C is now divided by 1+ pg,(A,p)/[uQ)—pl. We 
shall see in the next section that this factor tends to a constant limit as p approaches the singular 
point us 

If the ratio u,=y,/y is assumed to be known, yo can be calculated, and thereby y, becomes 
also known. The second equation (17) yields then y, directly, the third yields y,, and so on, 
that is, the whole system of equations unravels automatically. It was pointed out by Wick [1], 
and also by Waller [19], that the set of values thus obtained for the y,, diverges, in general, 
unless the correct value of u, was chosen at the start. Thus the solution of the system (17) 
takes the character of an eigenvalue problem. (This property is common to all systems of 
which the first n equations contain n+ 1 unknowns). 

A simple trial and error procedure on the choice of u,=y;/yYo in (17) is impractical in the 
critical range of \ where the directional distribution is sharply peaked forward and the successive 
coefficients ¥),Y.., . . . are nearly equal, so that the trend of the sequences y,; evolves exceed- 
ingly slowly. In other words, the angular distribution cannot be characterized in this range 
by just a few Legendre coefficients. Therefore Wick [1] replaced the sequence of the y,’s 
with a continuous function and the expression [(¢+1)y..,+/y..)/(2l+1)—y: in (17) with the 
second derivative of the function. In this manner the problem of determining the value of 
Us=yi/Yo was reduced to a problem of differential equations. An approximate determination 
of the eigenvalue by variational methods offers little difficulty and was carried out by Wick for 
a constant yw (A) and for a linear u(A) [1]. 

Whereas variational methods are well suited for deriving rapidly an approximate eigen- 
value, they tend to vield poorer approximations to the eigenfunctions and generally do not 
appear convenient as a basis for higher approximations. In fact, Wick effectively characterizes 
the angular distribution of radiation by a single adjustable parameter, in the first approxima- 





on. A substantially better approximation is required to take into account more general 
ariations of uw (A) as well as variations of the kernel & (\’, 4). Extensions of Wick’s procedure 
) this direction have been developed but they appear cumbersome. 

When the directional distribution is sharply peaked, as it is in the critical range, \~A,, it 
an be characterized conveniently by its set of moments, which in this case tend to converge 
ar more quickly than the set of Legendre coeffcients y,. The moments, mo=yo, m;=giYo 


rere * (1—u,)"y(p,u;,d)du are linear combinations of the Legendre coefficients. 


Jinx 
Determinations of g, by trial and error procedures based on calculations of the m,’s should be 
nore practical than procedures based on the y;,’s. 

Spencer eventually developed an efficient iteration procedure that permits to determine 
¢, by trial and error, utilizing only a very few—as few as two—of the equations for the 
moments [6]. This method has superseded all previously attempted approaches. 

One point, which emerged in the course of earlier attempts, might be worth mentioning 
because of its possible more general interest. The testing of trial values of g,(p,A) becomes 
particularly easy when the obliquity of penetration is less effective than the variations of 
u(A) in shaping the X-ray spectrum, that is, according to section 2.3, when d log y/dA\<1. 
In this event, an incorrect initial assumption regarding g,, not only causes the sequence of 
moments ™;,M, ...M, .. . eventually to diverge, but quickly to assume a grossly erratic 
trend. This critical instability appears to be a general property of systems that may be 
called “strongly nonself-adjoint.”’ It is discussed in appendix B. 


7. Effects of Small Deflections on Special Analytical Solutions 


As mentioned in section 6, Wick did the basic work on the effect of small deflections on 
penetration [1], for two typical situations, namely, constant u(A) and linear u(A Either of 
these assumptions about yu, together with the assumption (11) of a constant kernel /(\’,d), 
makes it possible to separate the direction variable u, from the energy variable \ in the transform 
equation (8). 

Specifically, the homogeneous part of the equation thus simplified has solutions that are 
products of a function of \ and of a function of (1—u,)p/|u()—p]. This to say that the diree- 
tional distribution maintains a constant shape as \ varies, but it contracts in width, as u(A) — p 
grows smaller, in proportion to [u(A)—p]/p. As a result, the factor 1+ pg,(A,p)/[u(\) —p] 
in (21 remains ind¢ pe ndent of dh and P. Therefore, the dependence of the photon 
number on \ and p, irrespective of direction, which is given by (21), has the same analytical 
form as the straight-ahead approximation (13b), except that C is replaced with the constant 
Ci\l+pa:/(u—p)}. There is a whole set of solutions of the homogeneous equation, with 
different values of PI/ (u— p 

Wick superposes solutions of this set so as to construct a solution of the inhomogeneous 
equation with the 6(A—d) source. The component solution with the smallest value of 
pqi/(u—p) predominates over the other components as p approaches 4, and as \ increases 
through the critical range \~,. Therefore, the singularity of the Laplace transform is char- 
acterized by the variation of u(A)—p in (21) (or in (13b)) as p approaches yu, and by the smallest 
value of pgi/(u—p). 

The application of this procedure to the X-ray problem is straightforward when the pri- 
mary radiation is the most penetrating one (Ap >Am, w= Ho) [4], since in the very proximity of 
\) the variations of u(A) may well be regarded as linear, u(A) ~uo+ go(A—%), and those of 
k(\’,.X) may be disregarded. 

The applicability is less clear when the primary radiation is not the most penetrating one 
ho Am), that is, when (A) goes through a minimum in the critical range. However, one 
may surmise that the problem at the minimum (A~),,) can be treated by assuming the same 
value of pg,/(u—p) as though u(d) were constant. This assumption, which has not been dis- 
cussed previously, implies that the singularity of the transform is given by the straight-ahead 
formula (14b) with C replaced with the value of C/{1+ pg,/(u—p)} corresponding to a constant x. 
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We present here first a qualitative argument and then a more detailed treatment which 


confirm the correctness of the surmise 

The trend of the solution of the transform equation in the critical range depends on the 
variations of the attenuation coefficient u(A) and on the effect of small deflections. We want 
to show that the deflection effects predominate over the variations of (A), in the limit as p 
approahces u,, so that the deflection effects may be treated as though y(A) had the constant 


value yu, 
Obliquity of penetration effectively impedes the build-up of photons if the value of 


u—pu, in (7) becomes substantially larger than the corresponding ‘“‘straight-ahead”’ value 
u—p. Therefore, the main build-up is confined to photons for which 1—u,S (u—p)/p 
Compton scattering of photons within this range of directions corresponds to wavelength 
shifts 6A<1—u,S (u—p)/p 

On the other hand, variations of the absorption coefficient that affect the asymptotic 
behavior must be of the order Au~u,,—p. These variations may be described, in the neighbor- 
hood in the minimum, as Au~}u, AX? (where 4,,=(d?u/dd*),,). Therefore, the wavelength 
change Ad required to bring about a substantial change of « is of the order of magnitude 
[2(4m—P)/Hml? 

Hence, in the limit for 4,,—p->0, this shift AX becomes proportionately much larger than 
the shift 5\ required to bring about substantial effects of obliquity. 

Proceeding beyond this qualitative argument, a complete treatment of the problem would 
consist of solving the transform equation (8) by a procedure of successive approximations, 
whose first step should V ield the Wick-type results. Efforts in this direction have yielded 
rather cumbersome developments. In practice such a complete treatment has been made 
unnecessary by the development of the semiasymptotic numerical method [6] (see section 2.4), 
which requires, as an analytical complement, only a knowledge of the singularity at p=u,, 
u,=1. One possible approach to the remaining problem would be to seek a special solution of 
the eq (8) for the case when the primary radiation is not the most penetrating one, that is, 
when Ay< Am, Hem, Since Wick has already solved the problem for \y>A,, uso. Whereas 
Wick assumed a linear variation of u, near \X=», we should assume here a parabolic variation 
near A=\,,, that is, u(A)~u,.+44,(A—A,,)?. Instead of following this approach, we shall take 
a somewhat more general one, namely, to seek special solutions of the equation valid “in the 
asymptotic limit” 


[u(A) p) Pp >), (22 


that is, for \~A,, up~1. This will be done without utilizing any initial approximation about 
u(A) or A(X’,A) but the approximations considered above will result automatically from the 
consistent application of the limit process. 

The first step of our treatment consists in choosing instead of u, a new direction variable 
that measures the effect of obliquity on penetration. Equation 21 shows that this effect is 
proportional to the mean value of p(1—w,)/[u(A)—p]. Accordingly, we replace the unit space 
vector u=(0,¢) with a homolographic projection variable that is represented in plane polar 
coordinates by the vector 


v= ((2p(1—u_)/(u—p)]'””,¢) ~ ([p/(u—p)]'"08,¢). (23) 


(The approximate equality corresponds to the small angle approximation 0< <1, which we 
need not utilize at this point.) The inverse transformation of (23) is 


$e*[u(d) — p)/p, u(v, \)=(arcos u,, ¢). (24) 


It is also convenient, though not essential, to renormalize the dependent variable y(p,u,,d) 
by multiplication with u(A)—>p, that is, to take 


w(v,r) =[u(A) —ply(p,uz,d). 
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[he dependence of w on p will not be indicated explicitly 
The transform equation (8) takes now the form 


AN KNX [ ve’s 1—u-u— dA+- Nd’ )wi(r’, rd’) + TL 6(A—A,) Fu 


ap [ ae’ A— A[X.0,07], A) wr’, A— A[A,0,0']) + 15(A— dy) FC — be? (y P), 26) 


where A stands for 1—u-u’. Since w is given by (24) as a function of v and \, and similarly 
u’ as a function of v’ and \’=A—A, A(\,v,v’) is defined as the root of the equation 


u(v.d)-u' (vd 


We are now interested in the solutions of (26) in the critical region X~A,, w(A) ~ w, and in 
the asymptotic limit [u(A)—p]/p~0. Following Wick, we seek here solutions of the homo- 
geneous equation corresponding to (26), that is, of 


(1+ 3e*)w(r,d) = (2p)! | dv’k(A— A, A) we’, A— A), 

and leave it for the next section to superpose solutions of (28) in such a manner as to fulfil 
the inhomogeneous equation (26) at the source where A= dp. 

Our knowledge of the straight-ahead solution and the discussion of section 6 on the effect 

( 


_— x 
of small deflection suggest that we split off from w(r, \) a factor exp of dd’ /[w(r’) Ply» 


where C remains to be determined, with the intention of showing that this factor contains the 
whole singularity of w in the asymptotic limit. Thus we set 


te fn On i 
wir,r x(v7,A exp yf dx [u(r } Pi 


— *” 
y(p,u,,) =[u(X) — p]'x(e@,A) exp | dd’ /[u(’) pit: 


Equation (28) becomes 


2aep)-' [deren A,X) exp) c| ? AN’ /[w(r’) pl, x(v’,r), (30) 
. A—A 
and we shall show that in the asymptotic limit A vanishes, so that #(A—A,d) equals C, and the 
exponential becomes independent of \. As a result x becomes, in this limit, a nonsingular 
function of v only, so that the exponential in (29) does indeed represent the singular behavior 
of w The values of C will be fixed by the condition that x (v) converge for large v. 

To carry out this program, we begin by studying the equation (27) which determines A. 
If we substitute u and u’ from (24), (27) becomes 


A 1 [1 = u(r) p e]f | lr A) pP | 
2 P 2 p 
fi u(A)—p p— i HO P 2 1 w(A—A) ef1-j =P y's my 31 
(2 p 4 p 2 p 4 p ) 


Terms of higher order in [u(\)—p]/p will be disregarded henceforth in the asymptotic limit, 
Such terms should be treated by a power expansion in case one intended to proceed to higher 


lll 








approximations. This procedure has the same effect upon (31) as the “small angle approxima 


tion” indicated in (23). Equation (31) reduces then to 


1 u(d)—p §— 2, w(A—A Pe |[ \) "] ‘2 en’: 
2 p i u(A)—p u(A)—p \ 


Since this equation indicates that A tends to vanish in proportion to [u(A)— p]/p, in the 


asymptotic limit, we replace A with the variable 
E=Ap/[u(d)— pl, (33 


which we surmise to remain finite in the asymptotic limit. The dependence of u(A—A) on A 
will also be indicated through a parameter which tends to a finite limit, namely 


n= (wr A—A)]/pA=[w() — w(A— Ea) — pl/p))/ Ele) — p] 
u(A—2 u(r np. 
With these substitutions, (32) reduces to the quadratic form 


20-v"*)nlé+[(e? +e" )?—40-0’ |=0. 36 


This equation depends on £ not only explicitly but also through 7. However, if the right side 


s 


of (34) is expanded into powers of [u(\) —p]/p and only the first term is retained, one finds that 


n du/dx Pp u(A)/p 


independently of & Therefore, € is given by (36) as an algebraic expression £(0-v’,4/p) which is 
finite and independent of [u(\)— p]/p in the asymptotic limit. Owing to (33), A vanishes then 
in the asymptotic limit, as we had surmised. 


Therefore, in this limit, A(A—A,A)=A(A,A)=C, the exponent of (30) becomes 
AX lu \’)—p] (C/a) In (1 — Eu P|. 
and the equation (30) reduces to 
C/2xp fac’ (1 E (0,0 ,u/p)u/pl”*x(e’,r 


The kernel of this integral equation depends on \ through g. Within the narrow critical 
range of photon energies, \~A,, in which we seek an analytical solution, w (A) can be expressed 
as a power series 


ep + e(A—A,) +1/2p.(A—A,)? +... (40 


Accordingly, one may seek the solution of (39) by expanding both the kernel [1 — &(v,v’, x/p) a/p\C/i 
and the unknown function x(v,A) into powers of (A—A,). Since the critical range around \ 
becomes infinitely narrow in the asymptotic limit, the singularity of y(p,u,,A) (at p= p,,u,=1) 
depends only on the zero-order term of the expansion into (A—A,)". Therefore, for the purpose 
of studying this singularity, we may set 


X(V,A) ~X(P,As) = Xo”), (41 


where xo(v) fulfils the eq (39) at A=A,, and with p=g,, that is, 


(Cc 2mu,) {dot £(0,0 .u My) us]°/"* xo(v’). 
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this manner we have reduced the problem to the solution of an equation in the single vari- 
that is, we have achieved the same separation of the energy and direction variables 
The only difference lies in having achieved the separation 


ie 
and r) as Was achieved by Wick 
a limiting procedure rather than by initial assumptions regarding the variation of uw(A) and 


‘constancy of k(A,A). Notice that (42) holds whether A, coincides with the primary wave- 
neth A» (that is, for X, A< \,,) or with the 
ivelength of minimum absorption A,, (that is for e<A,.: 
), the kernel of (42) reduces to 


\,,), With the longest wavelength \ (that is for dp: 
). In the last event im —O and, 


cording to (36) and (37 


S , 
» v 


s, -9 exp {[—(C/u,,)€(0,0’ ,0)] exp 


[1—E(0,D’ ,5/us) Wo! Ms|' 
Equation (42) is equivalent to the Wick eigenvalue equations and has finite solutions x9(v) 
nly for special values of c. As in the standard eigenvalue problems witha single independent 
ariable, the successive eigenfunctions x{”’(v),x6?(v), . . . xi (r) have 0.1...n.. 
odes. The corresponding eigenvalues cw , c a c de depend monotonically on 7. 
he eigenvalue C is the largest one, algebraically. Each of the eigenvalues 0 characterizes a 
ype of singularity of y(p,u,,A). As mentioned before, these singularities coincide with those of 
he straight-ahead solutions | 
~. The solution of the inhomogenous equation (8) or 
section, as a superposition of components, each of which is a solution of the homogeneous equa- 
tion, with different values of C. We anticipate here that the component with the largest value 


9’) is clearly the one with the sharpest singularity and therefore the one which predom- 


14) in section 5 except for the replacement of the constant C with 
26) will be constructed, in the next 


{ 


of © in 
inates and which alone matters in the asymptotic limit. 

The equivalence of (42) with the Wick equations is immediately apparent when A 
Um, sO, that is, when the kernel of (42) takes the form (43), so that 


. l — 
E +50), r IC’ 2x uml | de’ exp | = \C/me jijo-—o ‘h A 


This equation coincides with (W53) [1], as shown by substituting 


L MUk)28, C/lum=2k/M, Clum=(M+1)2/2M 15) 


») 


c 


When 4,0 no better procedure is seen, to show the eigenvalence with the Wick equations, than 
to backtrack from (42) to an equivalent homogeneous equation in the variables \ and u, and to 
handle this problem by Laplace and Fourier transforms in the manner of Wick. One may then 
identify xo(v) with a transform of solutions of Wick equations according to the formulae 


Ria = : ; 
Xo) (us Qrius) | dt expt -~ # f([2tis us)'/?v) 


- . 7 oe 
(u,/8x71) dt | dg exp (t—it'*q-v) t'~©"ed([u,/2u,)'"@) 


e ' 


(t,/S27) | dt Jig exp (t—it'g-0+ q?/u)t!-©*U ((us/2us]“Q). 


The functions f, ® and U obey the equations 


[e-i.— + u(d du) + ina? [f(w) (C/aw)f du’ exp ( u--u’ *) flu’ . 


}gmd(d da)? + a7'(d de)\—d tod da) +C exp ( - 07) (6) Co(e), 


' ‘ ‘ . re ' > 
gal (a do)*+a-'(d do)\+5 hes +C exp ( 0°) — 5 feo ey U(e)=CU (a): 
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All of these equations obtain when yu,=0 as well as when uv, #0, and for u,<0 (that is, for 
\.<A,,) as well as for u,>0. In particular, (47) coincides with (44), for u4,=0, if w is replace: 
with (20/u,,.)'20 

The self-adjoint differential equation (49) has the form of the Schroedinger equation fo: 
the motion of a particle oscillating radially in a evlindrically svmmetrical potential It IS 
identical with (W105) [1] as shown by performing the substitutions 


V/ 2M, mM Mi 24 V/, (’ vi ; p 


It is also identical with eq (V4) [4] as shown by replacing o? with o7/2 and C with —@Q/g,. 
Wick has estimated the largest eigenvalue of his Schroedinger-tvpe equations by a vari- 


ational method When xg 0 he finds 
51 


This eigenvalue reduces to C as u,,/C approaches zero. The ratio C/u,, represents the differ- 
ential cross section for Compton scattering, for \’ =A, expressed in units of the mmimum value 
of the total cross section. When 4,0 the largest approximate eigenvalue is 


where the parameter r indicates the solution of the algebraic equation 
l ae Meh / 


When the logarithmic derivative |g,/u,| is large, one may solve (49) by a perturbation 
method starting from the solutions of the quantum mechanical harmonic oscillator problem 
and treating e=yu,/u, as a small quantitv. The result, given in (V6) [4], is 


( 


with 
(1+4(C/i, el. 


In practice uw, gu, is seldom very small, and the range of application of the simple perturbation 
method is limited. The investigation of eigenvalues proceeds more effectively by the simple 
variational method of Wick or by the methods that were developed by Spencer [6] for the numer- 
ical solution of the semiasymptotic problem. Appendix C contains a description of these latter 


methods and tables of eigenvalues 


8. Effect of Obliquity of the Primary Radiation 


We resume here the discussion of section 2.6 regarding the influence of the source direction 
on the penetration law 

In the iast section, we have studied the singularity of special solutions of the transform 
equation at p=ys,. These special solutions obey a homogeneous transform equation shorn of 
any source term. It was understood that the solution y(p,u,,d) of the complete, inhomogeneous, 
transport equation (8), in the critical range \~, and in the asymptotic limit, should be con- 
structed as a superposition of special solutions (297) 


-_ "> 
y(p,u,r >, @y[u p| “ ver exp)" dyn’ [u(X’) -p 


*)r 
~Evelu(d) —pl-'xs” (0) exp} 0 | dN’ (u(x) —pl}- (56) 
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e x, (”) and C” represent, respectively, the n-th eigenfunction and eigenvalue of (42 
function x'(v,A) in the middle part of (56) represents the complete expansion 
Dx,” (vr) (A—,)” which begins with x," and can be carried to higher order of approxima- 

1 if necessary 
As p approaches y,, the factors exp) 0 (r’) pit become infinitely large, provided 


{ 


0. Wick pointed out that only one of the eigenvalues C, namely C®, may happen to 


positive. At any rate the term n=0 of the sum (56), with the algebraically largest eigen- 
uue C®, will increase fastest as p approaches u, and will predominate over all others in the 
.ymptotic limit. Therefore, the singularity of y(p,u,,s) depends primarily on the behavior of 


*A 
pyc dy’ [u(r’ pl, as anticipated in section 7. The behavior of this factor is described 


14) in section 5. 

In addition, one must consider the factor [u(\)—p]~', if w(A) ~u, and especially the depend- 
nee, if any, of the coefficient a) on p. The coefficients a, of the superposition (56) are deter- 
ined by the inhomogenity of the transform eqn (8), that is by the characteristics of the radia- 
on source Therefore, these characteristics may influence the singularity of YC pu d), and 
ence the penetration law at great depths, through the dependence of ay on p 

It has been pointed out in section 2.6 that the obliquity of the primary radiation has much 
ess influence on the penetration for high source energies, \y< <\,, than for Ay >A This im- 
vortant difference shows up clearly in the process of determining the coefficients a, of the ex- 

pansion (56). 

If the source energy is considerably larger than the energy of maximum penetration, one 

vill probably have to resort to numerical integration of the transform eq (8) from the source 
wavelength \y at least up to some wavelength A, somewhat lower than \,,. If A A,., and one 
s interested in the distribution of photons of energy higher than the energy of maximum pene- 
tration, then A, equals the wavelength of these photons and one must choose \;<\,. Whether 
AoA, OF Ay=A,, We define Ay so that Ay<CAySA, <A. A superposition of analytical solutions 
valid in che critical range \~A, should then be fitted to the numerical solution at A,. Under 
asvmptotie conditions, (u,— p)/p: 1, the fitting-wavelength A, can be so chosen that u(r, 

u u,—p, that is so that u(A)—p can be safely replaced with w(\)—x, for X<Ay,. When this 
is so, the photon distribution at A, varies no longer as p approaches yu, still further. Similarly 
{co “an [u(X’) pit, extended from A~ A, to A, no 


longer depend on p in the proximity of \,, where p may be effectively replaced with 4 


the analytical solutions x‘"(v,A) exp 


Under these conditions, the fitting of the superposition of analytical solutions to the actual 
photon distribution at A,, that is, the determination of the coefficients a, of (56), gives a resull 
independent of p. Any change of obliquity of the source radiation causes a change of the 
directional distribution of the photons at A, and, therefore, a change of the coefficients a,. The 
photon distribution through the critical range will vary with p in a manner which depends on 
the intensity ratio of the component solutions a,x,(r,A) exp of (56) and, therefore, on 
the obliquity of the source radiation 

However, the intensity of the lower energy radiation depends only on the n=0 term of (56), 
once p is sufficiently near to u,, and therefore it will vary with p only through exp 


— Te F . . , or 
( AN’ /[w(r’) P\y, independently of the ratios among the a,’s. That is, the source obliq- 


lity influences the intensity of the lower energy radiation through the same factor a) independ- 

ently of depth of penetration, provided this depth is adequate to insure fully asymptotic 
onditions. 

A quite different situation prevails when \)>A,, since the coefficients of the analytical 

olutions in (56) must be fitted to the directional distribution of the source radiation right at 

d,, in the middle of the critical range. Here, the dependence of the a,’s on p has been worked 

t by Wick for an isotropic source [1], and in reference |4] for a monodirectional source 
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General formulae equivalent to a determination of the a,’s will be given here below, but the 
are of very limited usefulness because of their complication and, especially, because they imp] 
a knowledge of the directional eigenfunctions xo"'(r). A practical solution of the effect o 
source obliquity under these conditions will probably require much additional work. 

An analytical solution for \)>),, equivalent to (56) and fitted to the directional distr 
X-ray source, can be obtained performing a Laplace transformation of the energy 


bution of the 
variable \, according to Wick [1]. This is a solution of the transport equation (8) schematized 


in the following manner 

(1) w(A) is taken as po+ Ho(A—A), 
(2) k(M A) is replaced with C, 

(3) u=(0,¢) is replaced with a “small-angle approximation” vector 8= (s,g), with O<s< @, 
such that wu, may be replaced with 1— 4s’, and u-u’ with 1—}/|s—s’?. For small angles, one 
The solution of the equation is then written as 


sets sn? 


. = 1 ¢ bo 
VPs =z: | dae | dhe! vinketiatian 


= DWE! | 
_ a : I ” dle “[at . ] D.(»[ 14 bol }) (; 
Ho—P Jo Mop Ho— Pp 
Here U’, and ©” are the eigenfunctions and eigenvalues of (49) and 
é kits DS F , . 
is, in essence, a coefficient of the expansion of the directional distribution of the source—in- 


dicated by f(u,) in (8)— in terms of the eigenfunctions of (47). Some simplification is attained 
by considering only the X-ray flux integrated over all directions, since in this case (57) yields 


; | r U7,(0) - jiol VO lio-1 ] 
d »8,A dne?™—o) S,, : | di “ft ] dD, a. . 
| 81/(p Sei. | n Fue—>) . ¢ a ( n — ) 


(59) 


Much further simplification obtains for two special types of source, namely, (a) mono- 
directional perpendicular to the source plane, f(1— }s*)=4(8), and (b) isotropic, or practically 
isotropic, f(1 $3?) ~f const., since (58) reduces, respectively, to 


D(a) = J kdke—*ok?/2P 17 (ke), (60a) 
2r Jo 


and 
D, (a) 4L,(0)f. (60b) 


Whether 7), is constant or proportional to its argument, the two integrals in (59) may be 
treated as a Laplace transform followed by its inverse, so that (59) reduces to 


° U0) « ” YE /in-1 6° Wei ts 
dsy(p,8,r) =>, : 1+ (A= Yo) F kdke—#o*?/2P [7 (he) (6la) 
4r dd wo—p Mo— P . ore 


bo ( A— No) oO /sig—1 > 
; (61b) 


f day(p,s,r) UO l =| 1+ bo — p 


respectively, that is, to algebraic functions of 4—p. (The p in the exponential of (61a) may be 
safely replaced with yo). 

For other types of directional distribution, the dependence of (57) or (59) on uo—p is quite 
complicated and even difficult to survey qualitatively, as had been anticipated. 
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9. Inversion of the Laplace Transform 


Following the discussion of the singularity of y(p,w,,A) at p= , in sections 5 and 7, we 
all now translate each type of singularity of the Laplace transform into a corresponding 
stribution-in-depth of the X-rays for large values'of the depth 2 within the medium. The 
anslation is done by means of the inverse transform formula (7 


. . 


l . l 
ru A) - | exp proy piu, r)dp eXp | Mel) >— exp (us p r| y(p,u A dp 62 
Ti i x 


2m, 


laving split off the factor exp(—4x,7) in the last expression, the remaining integral represents 
he buildup factor B(z,u,,) (to within a constant factor equal to the source intensity). Hence- 
orth, we shall be interested in the dependence of B on z for any given obliquity u, and wave- 
ngth \ and, therefore, we shall omit explicit reference to u, and X. Similarly, it will be con 
enient to indicate y(p,u,,d) simply as y(u,—p), so that (62) takes the form 


l 


") ~ 
2m, 


dp exp | (us P)rly (us p 3 


The preceding sections have called attention to two types of singularities of y(u,—p), 
namely: (a) an inverse power law, which we indicate as 


y (us— px (us—p)*" (64a) 


and (b) an exponential function of the inverse square root, which we indicate as 
y(us—p)x exp [Du,/(u;—p)|'”. (64b 


Equation (64a) describes the singularity when \)>\.,, and the directional distribution of the 
source is either perpendicular to the squrce plane or effectively isotropic, as shown, respectively, 
by (61a) and (61b). In addition, the same type of singularity arises for \y<\,, in the caleula- 
tion of the buildup factor for \<\,,, that is, when Ay<A,< A», in which case the straight-ahead 
result (14c) applies, except for the replacement of C with C®. Henceforth C® will be in- 
dicated briefly as ©. The value of K in (64a) is given by 

K=CU Ls | (65) 
for A,s< Am, according to (14c), and for A,>A,, with an isotropic source, according to (61b), 


whereas 


K=C/a,-44 (66) 


for \, >A, With a monodirectional source, owing to the derivative in (6la). Equation (64b 
deseribes the singularity when },=A,, and A>A,, which implies \<A,, in which case (14b 
applies, except for the replacement of C with C, and 


D=2rC? Um (d?y/dd*), (67) 


It may be added that the build-up of the radiation of maximum penetration A= ),,, for 
hich A,=A, also, involves an integration of the arctangent type of integral in (13b) over a 
inge of x/2 of the arctangent, instead of x. Moreover, in this case, the factor [u(A)—p]~' 
th) represents (u,—p)~'. As a result, (64b) must be replaced with 


Y(ue—p) (us—p)* exp [Fus/(us—p)}”, 


r(?/2y,,(d?u d»d*),.= D/A. 
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Notice that both K and D bave been defined as numerical coefficients whose value has t¢ 
be calculated for each medium and for the appropriate wavelength but is always of the orde: 


of 1. To emphasize the dimensional structure of A and D, we rewrite (65) and (67) in the 
following manner: 
a dy 
K= 7, k(a,»)(5- ) 


du /, 


Ore TCs). « 


The power-law singularity (64a), entered into the integral (63) vields the build-up factor 
proportional to «*, which was anticipated in (2). This result can be easily derived by replacing 
(u,—p)e with —t in the integral (63):* 


*i l 
dt exp (—?) {)-*-! 


l 
Birjx ‘ dp exp | eal er q-s&-!'—y, ‘ 
Sel. p exp [(us— p)t|(us— Py Oni. 
(See appendix C, tables 1, 2 for values of A, j,/us, C/ ms.) 
If desirable, one can also carry out the build-up factor integration by entering into (63 


not merely the singularity of y as given by (64a), but the expression of | dsy(p,s,d) as given by 


(Gla) or (61b). We may write 


ig io(A— Ao) 
Bir) 5. dp exp |(uo—p)2] (uo— p fi +* | , 
2m J —ia Mo— Pp 


and proceed, for example, by expanding the binomial 
r K ’ ‘ y 
' ol A — Ao) a rr + v) Mol A— Ao) 
n . . , 
Mo P —_ r(- K )p! Mo Pp 


and integrating term by term, which vields 


, h(—K+») : , , 
B(r)x >) =. Rt | My A—Ap)z]’ Fy K, l 9 MolA Aol Z) Lox - Mold Aol zr). 


Here ,/; indicates, as usual, the confluent hypergeometric function whose power expansion is 
| L_x« indicates the Laguerre function, to which the ,F; reduces 
The binomial expansion (71) is actually unnecessary 
Equation (72) 


represented by the =,, an 


when its second parameter equals 1. 
since (70) is directly recognizable as an integral representation of ,F, or L_x. 
was derived from a schematized treatment of the straight-ahead problem in (I4) [2]. 


The singularity (14b) for A,=’,, entered into the integral (63), yields a build-up factor 


represented by a transcendental function of unusual type, namely 


x 


| . 
Bir) =A oP exp { (um—p)t+[Dun/(um—p)|” 


—_ 
io 


wi ff . 3¢y 
Din(Diumt) ams ja exp { (¢+t~'”)(Du,z)'®} = Dun(Dunzt)"?G((Dunz)'”). | 


Here (= (um —p) (Du»)-**z'*, the contour A loops around the origin in the positive direction 


and 
_ l ; . 
G(z) ad Re exp [(¢+¢-'/*)z]. (74 


f the integral as a I'-function is found, for example, in reference [20], last formula on p, 245. 
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* The expression 





For large values of 2, the function @ can be evaluated by steepest descent method in the 


ximity of the saddle point, which lies at ¢ 


iis formula, entered into (73), vields 
B d oc (Du, (Dnt } 
his build-up factor is proportional to 


H=2-*3D) 


5(2 ar) 


® exp [/7(u,7)'"), as anticipated in (3 


2-*8 The result is 
“exp (2 

exp {2 3 (Dy md 

) with 


3(9°C*/2u mit m)' 


his result was derived in (II 13e) [3] for the straight-ahead approximation, that is with C in 


he place of C. 
The case A=A,=A,, 
itive of B(x) with repsect to J)! 


The function G(z) has been tabulated throughout the range 
ntegration (74) numerically along the steepest descent path Im(¢-+ 


(See appendix C, table 3 for values of u,,, 7, D, and (). 
with the singularity (64c), may be treated by taking the second deriv- 
*. since this operation generates (64c) from (64a), and then 
eplacing )) with the appropriate value (68) of F, namely ))/4, and H with /7//2° 


0.63 1 
0 by carrving out 
g-"*) = @ 


the 
However no 


application of this tabulation has been required because Spencer [6] found it convenient to 


evaluate integrals of the type 


ae; | CP(Ho— P)* exp { (Hm —P)t + [Dim/(um— py)” 


by expansion into powers of [Dum/(um 
manner of (70) and (72). 


sion Was moderately rapid despite the large value of the argument. 


(78) 


p)|'? followed by term-by-term integration in the 


The method proved practical because the convergence of the expan- 


(The number of terms in 


the expansion corresponds to the multiplicity of scattering [6)]). 
I thank L. V. Spencer for many helpful discussions and for contributing the appendix C. 
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11. Appendix A. Connection Between the 
X-Ray Distributions From a Plane Iso- 
tropic and From a Point Isotropic 
Source 


The X-ray distribution from a plane isotropic source is 
with f(u,)=1 The 


represented by the solution of eq (9), 


equation for a point isotropic monoenergetic source is a partic- 


ular case of (4) in which one takes the source term S(r,u,d) 
5(r)8(A—X). The solution of this equation will be indicated 
as F(r,u,,d), where u,~wu-r/r is the radial component of u 
Following up the discussion of plane-symmetrical geom- 
etries in section 1, we notice that the plane isotropic source 
5(r)8(A—Xo) may be regarded as a plane integral over a 
distribution of point isotropic sources located at ro=(0,yo,20) 


5(r)b(A — Xp f dyed 265 (7 — To)5(A— Ne (Al) 








Owing to the linearity of the transport eq (4), the distribution 
“(4,u,,) will be similarly expressed as an integral over the 

distribution F(r,u,,A Consider an isotropic point source at 
(O, yo, 2%) anda point Pat (7,0,0 The distance from the source 

to P is r=[z? + yo? oe]? and u, is equal to 

[r? + yo? eo}! We have 


dyod zok ((z 


If we take evlindrical coordinates (p.¥.r), with p jo" 
and y= arcos [| ) 270 42)! 2) becomes 


2) cos ¥ 


p'|2 


The distribution F(r,u,,A) may represented by 
Legendre polynomial expansion 


F 
and similarly 
Yiz Us 
"2s 
Since the J, d ¥ P lu *p cos yl 


simply 2¢P,(u,) P ir p’| A3) reduces to 


- 


{ pdpk Fo p* | 


| ” RAR FAR) Pla/R (AG 
Zz 


The Legendre coefficients with l=0 represent the total 
intensity distribution at each point, irrespective of direction 
of propagation. Since Py=1, the relationship (A6) simplifies 
for l=0 to 

Yu(z,)= | RaR PRD (A7 


Me: 


Fylaz,r xr OY (z,A)/ or (AS 


12. Appendix B. Note on the Treatment of 
Strongly Nonself-adjoint Systems 


Consider the system of equations 
Bl) 


In the limit p=0 this system has a trivially simple solution 
since the first equation contains the variable z, alone and all 
other equations can be solved chainwise. If p is sufficiently 
smali, the terms pzr,,, may be regarded as small perturbations 
with respect to the trival case p= 0. 

Our method of proceeding is this: If one assumes a trial 
value of the ratio z,/r,, the first equation determines the 
absolute value of z,; and z,. Thereafter the n-th equation 
determines the n+ 1-th variable, chainwise Owing to the 
smaliness of p, any error in the initial estimate of x,/r,; which 


causes a wrong estimate of 2} a,,.7,,, causes a much larger 
mel 


error in the estimate of z,,;, and so on. A trial and error 
method with successive estimates of x,/r, may thus proceed 
quite rapidly, because the amplification of the error in 
successive equations soon causes the trial solution to become 
grossly unrealistic, as judged by some suitable criterion. 

Experimental calculations were made with the homoge- 
neous system 


Ln = (Ln—-3 + In-2 + Tn a+1)/2, (B2) 





assuming +, —0 for negative n. This system may be 
as representing a game of chance of the parcheses 
which case 2 must be positive. Successive vari; 
linear functions of the trial value of x,/r»9 with alter 
rapidly increasing slope. Calculation up to n=8 s 
show that the solution becomes unacceptable uw 
~0.325 

\s a physical model for the considerations prese: 
one may think of a chain of pendulums with asym 
or “coaster’’, couplings, such that the motion of ¢ 
dulum influences the next pendulum on its right 
but the one on its left only weakly. If the right-to- 
pling were nil, then the motion of the whole set would 
determined proceeding from the left to the right 
the right-to-left coupling is weak, one can still deter 
whole motion easily by approximation methods proce 
in the same way from left to right. Any error in tl 
estimate of the motion of one pendulum implies a lar 
in the corresponding estimate for the pendulums fu 
the right 

We call this situation “‘strongly nonself-adjoint”’ becaus 
is a typical feature of the self-adjoint problem wit 
metric coupling that one must consider the whole set of 
dulums simultaneously Similarly, we would chara 
the trivial case of a completely one-way coupling as “ext; 
nonself-adjoint’’. 

In general terms one might state the situation as fol! 
Self-adjoint systems are represented by symmetrical 
rices. Extremely nonself-adjoint systems are represented 
matrices all of whose elements vanish on one side of 
diagonal; these systems have simple solutions. Systems 
matrices having small-valued elements on one side of 
diagonal should generally be amenable to easy approx 
treatment, especially if these elements are confined to a st 
near the diagonal. 


13. Appendix C.* Determination of Wick’ 
Eigenvalue by Semiasymptotic Itero. 
tion Procedures 


This method for determining eigenvalues can perhaps 
be illustrated by proceeding step by step through the solut 
of a specific example. Thus, suppose we want to deter 
the eigenvalue of eq (44) for the case of Pb, that is 
C/im=2.90. To begin with, we backtrack a little and def 
a variable u=(C/u,)'?». In terms of this variable 
becomes 


= > ' i si . 
[? mts U* [xo(u) C/um) | du’ 5—e'?|u— u’ |? xo(u ( 
é d 25 


Next, we multiply this equation by powers of (u*/2), integr 
over the u space, and thus write down formally mor 


equations of which the first two are 
(Cum) Ko + Ki=(C/um) Ko 
(C/um) Ki + K2= (C/um) (Ki + Ko) 


where A, 1/2) | du(u?/2)"xo(u). If we write 


insert numbers, these two equations become 
EKy+ K,—2.90K,y | 
EK, + Ky=2.90( K, + Ky) 


Now the right sides of these equations are the mon 
a smooth, positive, monotonically decreasing (as 
function, which we shall call R(u). We shall assume 
reasonable representation for this function is R(u) =} 
. . " 219) 
From (Cl) we then have xo(u)=[&+ w/2}-'be#"™ 
If we take moments of this quantity, we readily ol 
relation 

l 
A 


Ki/Ko=9= 8) 769) — |p 


*By L. V. Spencer 


now establish iteration 


, namely 


FEi(—z)). We 
(3) and (C4 


an 


1.4556 
1.4444 
0.40723 
0.49284 


rence between successive —Es is converging by a 
4 at each step. Consequently, &* should 
small fraction of a percent, equal to 1.404 

0.484C. This Wick’s 

0.499¢ (See eq (51 

ype of calculation can be extended easily by taking 
int more of the moment equations (C2). A calculation 
ise of four moment equations vields an approximate 
with that of the variational method 


be, 
This 
variational 


compares with 


le in agreement 
than | percent. 
teration procedure for the case of a monotonically 
g attenuation coefficient can up in the same 
Our original intention was to calculate the eigenvalues 


be set 


TABLE 1 Values of i,/u, and 


The middle equation 
while the 
of the 


n’th approximation 
in (C5) is simply a rewrite of the first of eq (C3), 
last equation in (C5) is obtained by fitting 
assumed R(u) to the right sides of eq (C3). 
The function f(z) can be calculated from standard 
tables. For interpolation between tabulated values we 
the relation f(z +46) ~f(z)[1+-6(1+1/r)]—6 
the calculation 
° 1.90, 8 
Thus 


where n refers to 


moments 


easily 


make a 


0.5, and s®é 


we guess that g(° ] 


0.95 


To begin 
This yields ¢ 
(0.95) = 0.5863 


Continuing, 


1.4864 g =1.4931 


1.4136 : 1.4069 


0.40219 8“ 0.40111 


0.48632 /f(pE 0.48481 
of tables 2 and 3 by this method using four moment equations 
and starting from a variational method value It turned out, 
however, that the calculations agreed to 
within about a percent or better with the variational calcu- 
lations in all parts of the tables 2 and 3 and thus we decided 
to present the numerical values of the variational calculation 
results Obtained from eq (51) and (52 

The agreement between variational and semiasymptotic 
calculations argues strongly for an accuracy of 1 percent or 
so in the eigenvalue determinations. However, in practical 
siuations, A and H depend also upon first and second deriva- 
tives, respectively, of u These are difficult to determine 
accurately. There may thus be quite sizeable errors in the 
tabulated values of /7 and D (table 3 table 1 


semlasvmptotic 


and of Me Ms 


for various materials and 


Cp, 


energies 


Source 


energy 


10 


x 








of the ecponent K=C/,, of the build-up factor for various combinations 


Taare 2. Values 
of the parameters u,/u, and C/u, 


Values of J./u. and C/m, for various materials and energies are given in table 1] 


pe = 0.02 


0. 0200 0. 1005 37 2.63 13.00 
o1gsy 1052 1. 328 5.49 
O1s46 73 OUSS ‘ S71 
01617 Os44 A ‘A7 
olsal 0687 SSN 

oles 52 0506 5 

0335 . 1696 


TABLE 3 Values of um, H =3[xC?/2yu,(d?u/dd*),)]', 
D=4H?*/27 and C for various materials 


C 
(em?/g per 
Material Compton 
wavelength 
unit) 
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